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INTHODXJCTION 


IC Sohwarzschild’s classical woik on absorption and difliision in 
the sun’s atmosplieie, and the continuation of these investiga- 
tioiis through A S. Eddington, J H Jeans, E A Milne and othois, 
have rendeied the theoiy of radiative eq[uihbriuin a definite 
chapter of mathematical astiophysics The piohleins connected 
with this theory aie of the following type. Gaseous inateiial 
filling a certain part of space is subject to given moident ladia- 
tion.* What is the distiibution of light and of tempeiaturo in the 
medium, when in ladiative equilibrium ^ In the outer layers of 
a star, wheie the cuivaturo can be neglected, i o Avhore the 
material can be consideied as stratified in parallel jilanes, those 
problems appear in thou simplest mathematical form, The 
astrophysicists mostly contented themselves with an approxi¬ 
mate solution of these jiroblems Owing to a certain inherent 
beauty, however, they aroused also the interest of the rigorous 
mathematician It is the purpose of this tract to attempt a 
coherent lepreseiitation ol all that has been achieved in the 
diiectioii of a iigorous tieatment of those siandaid problems. 

The rigorous solution of a physical problem naturally pio- 
supposes an exact formulation of the physical assumptions im¬ 
plied in the problem. Has then the rigoious solution more than 
merely mathomatioal signifloaiioo, when these assumptions prove 
right only with a limited degree of approximation? The answer 
to this question is not always in the negative. The Sohuster- 
Sohwarzsohild model of a piuoly scattering atmosphere, or the 
Schwaizsohild model of a purely absorbing gray atmosphere, 
both in radiative equilibrium, represent for instance typical 
standard models whioh play a similar r61o in the theory of the 
stellar atmospheres as the ^mtormodiato orbits’ in celestial 
meohamos. The Milne model, being oi special importanoe as a 

* More gonoially, on ilio bounclaiy the lodiation satisflos given oonditionfl. 
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Simple limit ease (infinite optical depth) of Schwarzsohild’s model, 
may well be compaied with Hill’s jieuodio orbit in the hinar 
theory, It is well known what the detailed study of an inter¬ 
mediate orbit means for the computation of the actual oiU‘ 
A more detailed and rigoious treatment of the above standard 
models of ladiative equilibrium should thus deserve more than 
meie mathematical interest. 

Besides the Imown lesults (theory of Milno’s standard model j, 
infinitely deep slab with no radiation incident on its suifaoo) tho 
reader will also find more oi less detailed discussions of otlu'i’ 
models. Milne’s model with ladiation incident on tho boundai'y 
(insolation of a planetary atmosphere, reflexion ofleot in oloHO 
binary stars), the Sohuster-Schwarzsohild model with an 
arbitrary law of scattering, narrower limits for the solution of 
the Schuster-Schwarzschild standard problem, Milne’s model of 
a planetary nebula, and other problems. 

From the mathematical point of view, the mam foatm'o of the 
above models is that they are governed by integral equations, 
which aie hnear in the standaid oases (pure absorption and gray 
material, pure scattering) Its important pioperty, namely tho 
positivity of the kernel, has already been recognized by Schwarz- 
schild as fundamental for the rigorous discussion. It may bo 
mentioned that the reading of the book requires no special know¬ 
ledge of the theory of integral equations, since the positivity of 
the kernel allows all pioblems to be tieated m an elementary 
way Only Chapter iv presupposes some knowledge of tlio 
Fourier integral 

It IS a pleasure to the author to thank Prof Milne for many 
helpful remarks and the Syndics of tho University Press for 
accepting tins book as a Tract 

13. 11. 

Nov 1, 1933 



CHAPTER I 

FUNDAMENTAL PRINCIPLES AND PROBLEMS 

§ 1 TRANSMISSION OF RADIATION THRODGH 
ABSORBING AND SCATTERING MATERIAL 

The outer layers of a star are constantly exposed to an enoimous 
flow of radiant energy coining from the deep interior In working 
itself through the outer layeis, a part of the radiation is absorb od 
and Bcatteied. The scattered radiation will immediately be re¬ 
distributed in the diflteient dnections, the fieqitency being 
unchanged On the other hand, a paiticle heated up by the 
absorption of radiation re-emits temperature radiation of all 
frequencies 

As a fiist approximation, the state of matter and radiation 
may he considered stationary Therefoie wo must first find 
the conditions for a steady state The intensity of radiation 
changes in a definite way along the ray A part wiU be lost through 
absoi ption and scattering On the other hand, the inien sity gams 
again since the partioles on the ray emit scattered and tempera¬ 
ture radiation in the direotion of the ray This law of transfer holds 
separately for each froquenoy, 

Comhmed witli the condition for a steady state the oneigy prin¬ 
ciple yields another equation The net loss of radiant energy 
within an element of volume (I’adiation of all frequoiioies) equals 
the net gain of heat energy through conduction and convection 
plus the radiant energy liberated from other (sub-atomic) sources. 
Energy of the latter land must play a r61e in the interior of the 
stars m order to maintain their radiation over cosmioal ranges of 
time. 

When energy is transferred by radiation only, the net flow of 
radiation through each volume element vanishes. I’his state of 
affairs is usually called radiative equilibrium (in the strict sense). 
Modern astrophysics claims the outer layers of a star to bo 
approximately in radiative oqmlibnum. In the stellar interior 
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2 IPUITD A MENTAL PRINCIPLES AND PROBLEMS 

the energy hberated is no more negligible tlioiigh being small U i 
coinpanson to the ladiant energy itself In spite of this fact i t- 
has become cnstomaiy to speak of I’adiative eq^nihbrium (in 
wider sense) in this ease 

For the determination of the temperature distribution, oiit^ 
must Imow the laws relating the temperature of matter to tlit^ 
radiation When matter and radiation are in an onoloHUii>]i 
ICirohhoft’s and Planless laws hold. In the stars, howevor, th<' 
temperature vanes from layer to layer, so that those laws oaiuui tr 
be strictly vahd In the interior, however, the radiation is prae- 
tically enclosed, and even in the outer layers those laws luwo 
proved to be a very good approximation. 

In the steady state all C[uantities desoribing the state aro in- 
dependent of time The i adiation field is described by the inlansiit/ 
of-iodtation / (P ?) 

as a function of the point P, the direction ? and the frequency »».* 
Its physical meaning is the following. Let dcr be an element oT 
surface about P with the (sensed) normal n, and let day bo ail 
infinitesimal bundle of directions oontaming ? (surface elomonb 
on the direction spheie) If we write 

dq—\coB{n,7)\d<j .(1) 

for the projection of da upon a piano perpendicular to r, thct 
expression I,{P,,)dqdmdv .... (2) 

represents the radiant energy of the spectral interval dv which, 
per umt time, flows through da and spreads out in the sohd angle 
do>. Here is meant, the bundle of all rays through all points of da 
and m all directions contained in day dq is evidently the porpon- 
dioular cross-section at P of the bundle It is oonvoniont to havo 
a special name for the quantity 

I^{PyT)dqdv (3) 

We call it the strength at P of the parallel bundle of the spookai 
interval dv, of direction r and with the cross-section dq. 

We set dadv'^y^^^^dadv\Iy (P, r )oos (w, r)day, .(4) 

* The degree of polarization :s not analysed hoio, 
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where the integral is extended through all directions and whore 
n IS a given direction Let tZo- be an oiientated surface element 
about P and with n as its positive normal, The q[uantity 
represents then the net flux of radiant energy (of the spectral 
interval dv) through da and per unit time Energy flowing thro ugh 
da from its negative towards its positive side (cos > 0) is here 
reckoned positive, negative, however, in the opposite chroctions 
r(oos < 0) Let a, ?/, z be three orthogonal diiections, and let us 
write a, y for the three direction cosines of the normal n with 
respect to a;, y, z The equation 

cos {n, r) = o'oos (?, a)+ ^cos (?,?/) + ycos {'i ,z) 
shows that we have 


= cos ( 71 , x) g:„,a,+ OOS V) + {n, z) 

r e that g'v.n. component in the direction 71 of a veotor . 
This vector (P) is briefly called the net Jinx at P of the 
v-rachation. 

In a veotor field we have the Gauss integral identity 

= .(^) 

where F is a volume, S its surface, and n the outwaid normal of 
the surface element da 

For the total radiation of all frequencies we sot 


Z = 


■* CO 

0 


lydV, 



.( 0 ) 


We now come to the quantities that desonbo tho internotion 
between matter and radiation Let rj^, (P) bo tho mass coGjfictent 
of enmaton at P and within dv The quantity 

r)^,{P)d7ndo)dv .... .( 7 ) 

represents then the radiation omitted, per unit time, by tlio mass 
elomont dm at P, within dv and within tho solid angle doj In 
analogy to the above wo find it convotnont to call 

y]j,{P)dmdv .(8) 
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the stiength of a parallel bundle emitted by dm witlim dv l<^)i 
the total emission of all frequencies we set 

■n(P)^^\,(P)dv ... (fl) 

Let, furthermoie, (P) denote the mass coefficient of absorption 
for the v-radiation at P Along a short path fZs, the amount 
pv.ylyds of the intensity is absorbed, p being the density of 
matter at P Outside of astrophysics it is raoio customary to 
introduce the linear coefficient of absorption, pcf^, 

The mass coefficient ofscattei ing Oy (P) is defined in an analogoiiK 
way Soattermg weakens the intensity ly by the amount pOy lydn 
along ds pay is the Imear coefficient of soatteimg The radiation 
lost by scattering at P -will immediately be redistributed among 
the different direction issuing fxom P This distribution m 
described by a law of scatter mg 

.(10) 


where i denotes the direction of the incident ray Wo have, (jf 
course, . ^ 

Jyv(-P>n»')^^w' = 47r .... (U) 

The law of scattering is supposed to have the reciprocity projiorCy 

YviP>r>i')==yy(P,i',r) .(Ip) 

= .. .(U") 

wheie ~r denotes the chieotion opposite to r, lo the sanio 
amount of the soatteied radiation is sent into r and into — /, 
These conditions aie fulfilled under very general assumptions 
about the material The simplest case y, = l (unifoim scattering) 
has hitherto found the chief attention of the astrophysicists. It i«, 
however, necessary to consider inoie geneial laws too, for instance 
Rayleigh’s law , , 

yv-y(r,? )~|{H-oo82(^,?')} 

It IS seen to fulfil (IP) and (11’’), as docs also every law for which 
y„ IS a function of j cos (r, P)' 


All quantities (except the net flux) are non-negative, 
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§2 ABSORPTION AND SCATTERING OF THE MASS 
ELEMENT THE EQUATION OF TRANSFER 

We now compute the v-radiation absorbed and scattered witliin 
a mass element dm at a point P For this purpose let us hist con¬ 
sider the parallel bundle of all rays through dm which have a 
given direction u It is convement to imagine dm subdivided into 
many tliin columns of direction r Let ds be the length, dq the 
cross-section of such a column. The strength (P, ?) dqdv of the 
parallel bundle of direction f through that column, i e. through 
dq, is, in consequence of absorption witlnn the column, weakened 
by the amount pcf^I^dqdvds Since pdqds is the mass of the 
column, we infer by summation over all columns that 

c/yl^,dmdv 

is the amount that the strength of the bundle considered above 
loses on account of absorption within dm By strength of that 
bundle is meant hei e the sum of the strengths of all the above 
partial bundles, once taken when entering dm^ the other time 
taken when leaving dm (x.f,Iydmdv is the difforonoe of these two 
sums Now, eneigy is obtained by directional integration ol the 
strength Integration over all directions ? of the above loss of 
strength shows thus that the whole d!v-radiation travelling 
through dm per unit time loses the amount 

dmdva^{P)iI^{P,7)dw , ..,(12) 


in consequence of absorption 

In the same way it is seen that the amount 

d/mdvcr^,{P)^I^,{P,r)d^lO .( 18 ) 

is scattered within dm. We now show how this scattered radiation 
is redistributed among the difioient directions i'. Aocordrngjjo 
the above considerations, the radiation llowmg through dm per 
unit time and spreading out within a solid angle du) (containing r) 


loses the amount 


M — dmdv Qy (P) ly (P, r) 


through scattering within dm. Tins soattored amount B wiU bo 
distributed over the other diieotions r' aoooiding to tho law of 
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scattering. The part falling witlun a cone dco* about tlio dirootioii 


r' 18 thexefoie 


^y,(P,ry)dm' 


Summation over all cones c?aj, i.e integration through all inoidon t» 
directions ?, gives, tlierefoie, 

dmd<o'dvjl,(P,i)YAP.i,i")dui ....(14) 

as the part of the scattered radiation (13) which is lodistributed 
within the cone dco' If we now integrate this over all dirootioiiH 
we obtain, according to (11), again (13) as it should bo. 

We note for later purposes that the parallel bundle of scattered 
radiation emitted by dm into the direction r has the strength (tlicx 
notation for incident and scattered ray is interohangod) 

dmdv f I, (P, r') y, (P, »•',)) rfw'.(14') 

Mquakon of transfer Let us consider a light ray of direction r 
through P, and let P' be a nearby point lying in the direction r 
from P, ds=PP' We construct a thin cylinder about da with the 
two bases at P and P', and witli the cross-section dq, The strength 

lAP'ydqdv, .(15) 

of the parallel bundle of dnection r through dq consists thoii of 
thi*ee different pai ts Pnstly, the strengtli at P of the same Inindlo 
weakened by absoiption and scattering in the cylinder, 

{1-/3 {cfy + of)ds) (P, f)dqdvi 

secondly, the strength (8) of the iiarallol bundle omitted by the 
cyhnder in direction r, with 


dm^pdqds^ 

thirdly, the strength (14') of the bundle of scattered radiation 
sent by the cyhnder into the direction r On dividi ng the oqu ath n i 
obtained in this way through pdqdadv, wo got Sohwarzsoluld’H 
fundamental equation of transfer 
1 


P + Vv) da 


'^,(P,r)~/„(P,r), .(10) 
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together with 

(»„+<7,) (P,))=(P)+^ J/„ (P, /) y„ (P. >»•) do-' 

. .. ( 10 ') 

The derivative on the left-hand side of (16) is the diiectional 
derivative with lespect to P, taken in the direction r of the ray. 
On introduoing three orthogonal axes x, y^z, P~ (jf, y, z), and on 

etting (X = cos(r,aj), j 3 = cos(?‘, 2 /), 7 = 003 ( 1 , 2 :) .,.(16") 

for the direction oosines of ?, we obtain 


di di dJ 

ds ^dx^^dz 


.(16'") 


The quantity mtioduced by (16') may, aocoiding to 

Schwaizsohild, be called the Bigiebighe/it It is an aveiago of the 
two quantities , « 


having thus the dimension of an intensity. Tho Ergiobigkoit 
depends, m general, upon the diioction j. In the case of uniform 
soatteimg (y,,= l) or of isotropio radiation (tho mtonaily is in- 
dejiondent of r), however, it is a function of P only 

The equation of transfer is of great generality It does not 
require that radiation is the only mode of energy transfer, it just 
refers to tho part of energy appearing in the form of radiation 


§3 RADIATIVE JSQUTLTBBTUM. 
RADIATION AND TEMPERATURE 

If ^iTT€{P)dm denote tho heat eneigy hbenUed in dm, 1 . 0 , tho sum, 
not gam within dm of heat energy due to convection and conduc¬ 
tion pins tho radiant energy (of all fioquoncios) coming from sub¬ 
atomic souices within dm pei unit timo, tho oonsorvation of 
energy is expressed by the equation 

f ^'.ido = 4Tr f pedv, .(17) 

js Jr 

where V is an arbitrary volume, 8 its surface and n tho outward 
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normal on S. Aocording to (5), intogj'aiod through tho wholo 
speotrum, this is eguivaloiit to tho equation 

chv^ = 'J7rp6 .. ,.(17') 

It 18 of impoitance to imli/o that (17) or (17') ioIoir only to l!io 
total radiation of all fiequoncios iSoattoiing luiR no in/luouro 
upon the wave-longth, tho absoihod ladiation is, howovoi, 
re-emitted within other parts of tho Rpoctiuni, iiui)lying, in 
general, an exchange of energy botwoon tlio (lifToiont jiartH of I ho 
spectrum 

Integration through all directions r of tho o(|uation of tiannlor 
(after multiplication with p {</,, H- ct,,)) yields, according to (11), tho 

equation div^'j, = 4.TTpY}y-p(/^JIydto, ... .(18) 

the scattered radiation having droiqiod out as it should do. On 
integrating (18) through tho whole spoctruin wo olitain, according 
to (17'), the fundamental relation 

t; (P)=i J cho r (^') I, {P,r)du + c{P). ...(Ill) 

If (19) is multiphed by L-ndm, tho left-hand side liocoinos, ac¬ 
cording to (8), the total omission of dm, whilst tlio iirst term on 
the right becomes the part of tho total radiation that is absorhod 
by dm (19) means thorefoio tliat tlio radiation omitted ininus tlio 
radiation absorbed equals tho energy liberated. 

Since (18) is a consequonoe of tho equation of translor only, 
(17') follows, conversely, from (10) (17') and (19) are tlioiofoixi 
equivalent expressions of the ooiisorvaiion of energy. 

In the ease e=0, used as a first apjiroxiination in tho outer 
layeis of a star, wo speak ol strict ladialive equiUbtum. JDnorgy 
18 transforied by radiation only, It should bo notod that c ooukl 
bo negative as well as positive m tho gonoral oaso. WJiou, for 
instance, tho element loses heat by oonvootion or oonduotioii, 
subatomic energy being missing, e is nogativo.* 

* (17'), (19) possess, of course, general validity, as they OKpi'cas tho onoig)’ 
principle They are, liowovoi, of impoitanco only wlion ladiatloii is tho pilnnlpnl 
agent of onoigy tiansfoi. It is m this sonso that (10) is apokon of as tho equation 
of ladiafcivo equilibrium. 
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Pwe abaoiphon Qiay rmU’nal We consider the case o-^sO. 
Furthermoie, let ns suppose that the absorption coefficient is in¬ 
dependent of the wave-length, «,,==« (gray mateiial). This case 
IS of special impoi tance because of its simiihoity The fundamental 
equations (1C), (lO') and (19) become, after nitogration thiough 
the spectrum. 


1 

vw) 


di{p,)) 

da 


=j(P)-/(n 


( 20 ) 




(200 


Monoohiomakc ladzaUve aqtnhbiviim In the case of pine scat¬ 
tering, wo know that the wave-length remains nn- 

olianged Tlie scattered radiation is fully le-emitted with the 
same frequency and distnliuled through the dilTerent cliieotions 
Tins state of altans is often called monochiomatic radiative 
equilibrium The equations of transfer (10), (10') take, here, the 
form , IT (P \ 

))<««' .... ( 21 ') 

In this case, we have, of course, according to (18), 

div^',^0, . ..(22-) 

It 18 xihysioally plain that the case of purely absorbing gray 
material in stiiot radiative equilibrium is formally equivalent to 
the case of monoohromatic radiative oquibbriuin and uniform 
scattering, y„sl. The absorbed radiation is uniformly le- 
emitted in direction, m the same way as the scaitored radiation 
is umformly rodistiibuted m direction The equations (20), (20'), 
es=0, have accordingly the same form as (21), (21'), y,,== 1. Tliis 
foimal equivalence is of importance, since it allows us to doveloj) 
the same mathematical theory for the two physically different 
oases, , 

Local tli&nnodynamical eqmUlyiium, In the case of pure 
scattering, radiation has no relation to the tomporaturo. As soon, 
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however, as absoiption and emission play a role, we should have 
information how the radiation is related to the temperatuio of 
the matter 

Let = be the intensity of the v*radiation of a black 
body at temperatuie T, 

QhvlLT^l> 


wlieie h sigmfies Planck’s constant, h Boltzmann’s constant and 
c the velocity of light Poi the total radiation we have 

B=^r BJv^-T\ .. .(24) 

Jo w 


<7 being Stefan’s constant. 

When the radiation is enclosed between black walls, it obeys 


Kirchhofif’s law 


f)ll OCj, By 


(26) 


The condition of being enclosed is, of course, not rigoiously 
fulfilled in the stars, since there is always a radial net fiux of 
radiation Eddington’s perfect gas stax, however, has such a high 
opacity that the radiation is practically enclosed, the degree of 
accuracy being higher than in laboratory exporiinents Astro¬ 
physics uses the term ‘local thermodynamical equilibrium’ in 
order to express that the radiation behaves like an enclosed one, 
i.e that Kirchhofi’s law holds Even in the outer layers of a star 
this state of affairs has proved to be a useful appioximation 


§4 THE MAIN PROBLEM 

In the physics of the outside of a star, tlio following pioblem plays 
an essential role. The rachation field, in particular the radiation 
emergent fiom the suiface and its distribution in dueotion and 
wave-length, and the temperature distribution are to bo deter¬ 
mined when the coefficients of absorption and of scattering as 
well os the law of scattering and e are given 
This problem naturally represents only a jiart of the general 
problem of the determination of the whole physical state. Tho 
quantities regarded as given in the above problem aro not stxictly 
kno’wn in the stars Actually they enter other more or less known 
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physical laws {containing of course the conditions of mechamoal 
equilibrium) which, together ivith the laws of radiation, should 
enable us to deteimine the total physical state The excessive 
difficulties, however, connected with this pioblem as it stands 
compel us to select partial problems, of which the one foimulated 
above is one of the most impoiiant. Compaiison with observa¬ 
tional data of the solution of this problem, taken for diffeient 
choices of the given quantities, has led to valuable insight into 
the stiuoture of stellar atmospheres The following compara¬ 
tively simple special cases have hitherto found the chief attention 
of astrophysicists Firstly, the Schuster-Schwaizschild model 
of a puiely scattering atmosphoio Schwai zschild was the fii&t 
to approach this model with rigorous mathematical methods 
Secondly, SchwarzsohikVs (formally equivalent) model of a 
purely absorbing gray atmosphere (6 = 0) The important limit 
case of infinite optical doptli has (by apjiroximate methods) 
been oxtonsivoly treated by Milne, who also studied the spectral 
distribution of the emergent light The explanation of the law of 
darkening on the sun’s disk is one of the main successes due to the 
study of those models of radiative equilibrium, SchwarzsohikVs 
investigations have, furthermore, led to the result that the origin 
of the Fraunhofer lines m the solar speotiuin is due to scattering 
rather than absorption. 

So far as concerns the general solution of the above problem, 
the following remarks aio of importance. The equation of transfer 
(16) can bo considered as a first order diftereiitial equation for the 
intensity. On integiating it along a ray, with resjieot to the 
boundary oonditiona, wo obtain the intensity exj)i eased in terms 
of the Ergiebigkoit If tins is insei ted into (10'), a linear integral 
equation (duo to Schwarzschild) is obtained foi the detoimina- 
tion of tho Ergiobigkoit In the case of pure soattei mg, ^^ v “ ^» 
this equation completely determines the Ergiobigkoit In the 
photospheric layers of a star, however, absorption piocesses aro 
predominant, and tho unknown omission = onteis The 
solution of SchwarzsohikVs intogial equation detoi mines, in this 
case, the Ergiebigkoit, and therefore the intensity as a ‘funo- 



A-ft jj'TTNDAMENTAI/ PRINOIPLiaS AND PROBLEMS 

tional’ of £f^, for every ficquonoy j/,* Since By=^B^{B) is a 
definite function of the intensity appears thus a definite 
funotional of. B~B (P), z.e. of the temperature distribution. If 
the intensity thus cToierinined is inserted in the equation of radia¬ 
tive oquilibiiuin (I9)> the iight-hand side becomes a definite 
functional of P (P). The left-hand side 

«.,(P)B,{B(P)}dv, 

however, loprcsents a definite function of the iioint P and of 
P « P (P) Wo thus arrive at a dofimto functional—^moie pre- 
oisoly, at a non-linear integral equation for the deierniiiiation of 
P (P), i 0 of the distribution of temperature Once B (P) is 
known, we aie able to find and theroforo for every fiequenoy 
It is diiHoult to solve the integral equation in its full gonerahty, 
and the limitation to simpler spooial oases seems necossaiy. If 
the mateiial is gray, i e. if is independent of v, the equation 
of radiative equilibrium simplifies oonsidorably If, however, 
vanes with v, the mtegiol equation still remains non-lmear. In 
gray material without soattoring wo have the sinnilest case, since 
the integral equation then becomes linear, and the determination 
of B (P) aiipears entirely separated from the determination of 
the spectral distribution, the latter being a mattoi of direct 
integration. It should be noted, that the mathematically more 
general case, where and o-y are both independent of r, loads to 
the same formal integral equation (hnear) as m ease cr,, = 0. But 
the spectral distribution of light must then be determined by the 
solution of SohwarzsohikVs integral equation. 

The greater iiart of this tract will he devoted to tho mathe¬ 
matical problems arising from the two particular oases, pure 
soattering, and pine absorption in gray matoi'ial The first case 
plays, as mentioned befoio, a oluef rdle in the outermost layers 
of a star, whilst the second case has found chief application to the 
photospherio layers. 

* This deteijiunfttion of Jy ja ijossihlo by moans of tlio Neumann sonos for 
the solution of the mtegial oquation, Convoigonco and iimquoness ol! that 
solution are tieated m f 34. 
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So far as concerns tlie boundary conditions, the radiation in¬ 
cident at the snifaoe of an isolated star is zero, while the radiation 
at gieat depth becomes isotropic In the case of a close binary 
stai, the radiation of the othei component must bo taken account 
of Or, m the case of a ‘ planetary atmospheie * in radiative oquili- 
brium, the solar radiation must be considered. In all these cases 
the radiation incident at the surface is given Another type of 
boundary condition ooouis in IVIilne’s spherical model of a 
‘ planetary nebula’ in radiative equihbrmm. The inner face of the 
shell receives not only the radiation of the central star but also 
the hght coming from other parts of the inner face The boundary 
condition appears here as a relation between incident and 
emergent radiation. 


§6 SOIIWAliZSOHILD’S MODELS MATERIAL 
STRATIFIED m PARALLEL PLANES 


In the atmospheies of celestial bodies the material may bo con¬ 
sidered stratified in concentric spheres The radiation field having 
spherical symmetry too, the intensity I (P, 7*) becomes a function 
of the distance a from the centie and of the angle 0 between the 
radius vector and the direction ? only, I^I{a,0),0^0S'n' In this 


case we obtain 


dl M smOdJ 
ds da a dO 


When, instead of a, 0, the variables 

^ = aoo80, 'q^aBmO 
are introduced, wo have more simply 


ds 


Neglection of outvalute, Let us consider a layer a^^a'^aQ of 
the star. If, the depth — Wo being kept fixed, the radius tends 
to infinity, we formally obtain the limit case of a "plane slab 
Another way of getting this limit case is the following. We sot 


= ^0 » P + ^v) “ 


a’o being a fixed quantity, while J{x) is a fixed function of x, 
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0<x<Xq. On proceeding to the limit aj, i o. on making the 
shell thinner and thinner and making p (a,, + a^) laiger and larger 
m the same ratio (at homologous points), we obtain, as seen from 
(16) and fiom the above expiession of dljds, again the limit ease 
of a plane slab. This time the ladius has been kept constant. It is 
customary to neglect the curvature in tlio oiitei layers of a star, 
the mateiial being thus stratified in parallel planes Let x gener¬ 
ally denote the depth of the point P below a fixed layer, a being 
reckoned negative for points above that plane The dnection ? is 
oharaoteiized by two angles 6, ^ (we letain the second angle <f} in 
view of later applications), 6 being the angle between r and the 
direction of negative (outward normal of the slab), 0 g 0 g ir 
^ 18 the azimuth, i e the angle between the plane, containing r 
and the a-axis, and a fixed plane through that axis, 0 ^ ^ g 2Tr. 
We shall only consider the case whore the radiation field is the 
same in all points with the same x, 1^ becoming then a function 
of (a, 1 ) alone, We have 

d8=~QQQddx, .(26) 


thus yielding 



-COS&;r-* 

ox 


Tlie ScJma)ZscMld-Milne model. This is the classical case of 
puiely absorbing and gray material in local thermodynamical 
equilibrium. It is convenient to intioduce, instead of -r, the 
optical depth p-e 

T= I apdx 

J -oo 

below the suiface In the case of strict radiative equilibrium, the 
fundamental equations (20), (20') take then, according to (26) 
and (26), the remarkably simple form 

cos<??^^|^ = /(r,0)-J(r) .. ..(27) 

and J (r) = i J Z (t, r) dio^ ... (28) 

where, in the case of local thermodynamical eqmlibiium, 


rj 


'00 
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If oonveotion, conduction or subatomic energy souices are not 
negligible, we must add the teim e/a on the iiglii-hand side of (28). 
An important fact is that the absorption coefficient appears only 
indirectly, namely tbrougli the mediation of the new independent 
variable r It is to be remaiked that (27), (28) determine the total 
radiation field, in particular the diiection distribution of the 
emeigent light, without any hypothesis about temperatuie and 
its relation to the radiation We need it, however, as soon as the 
spectial distribution, or the temperature distribution, is wanted 
Let us introduce the notation 

f =f • f 

J I J 0<nl2 J - J 0>7r/2 

foi the integrals over the hemispheres of all outward and all 
inward directions We set 



7(t,/) oo8 0fZw, 



/(t,?) I coa01 dcD, 


The equations 
following from 


I, 


cos Odco = 


cos 0 I dw = TT, 


dco=mn0d0d{l}, 


..,.(29) 


show that F ^ and F_ lopresent averages of the intensity within 
the respective hemispheres, l^uithermoro, we set 

TTF = 7rF^ — 'n-F_==JI QosOdo) . ..(30) 


The flux vector is, of course, normal to the layers, and we hav6, 
according to (17), 6 = 0, 

^ = 7rF = const. . (31) 


Another important relation, due to Eddington, is obtained if (27) 
IS multiplied by cos 0 and integrated through all directions 
According to (30) and (31) wo got 

As- f Jcos^0da) = 7V + const. . ,..(32) 


The physical moaning of K is that it equals c/tt times the radiation 
pressure normal to the layer. 
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In the absence of strict radiative equilibrium, 

the fundu' 

mental equations aie 3 / 

cos 05 - = 7 —7, 
dr 

. ..(33) 

J=^{ldco + -. 

It? J a 

(34) 

Instead of (31) we have ^ 

dr oc’ 

... . (36) 

and (32) must be replaced by 


II 

... .(36) 


Another simphfication of the model is effected by assuming the 
slab to have infimte optical depth, 

O^r <00. 

This ooriesponds well to the conditions •vvitliin a star (largo 
opacity in the interior) We formulate then the first standard 
problem, treated approximately by Milne: 

Problem I J (t, ?) ^ 0 and ^ (t) ^ 0, t < oo aie to be deter- 
mmed from (27) and (28), the radiation incident at the surface, 

7(0,?) ^0, 6 > 77 / 2 , 
being given, as well os the flux constant F. 

The incident radiation is zero for ordinary stars; positive, 
however, for close binaries (or in Milne’s model of a planetary 
atmosphere subject to insolation) For an ordinary star, the flux 
constant is derived from the observed effective temperature Tg, 
At the surface we have F = F ^, this being an average of the in¬ 
tensity of the emergent radiation. If this were due to a black 
body, the intensity would be ctT'^/tt, T being the sui’faoe tem¬ 
perature In general, 

.(37) 

defines the effective temperature Tg , 

The relations (27) and (28) being hnear in 7, t/=a J5, the solution 
of Problem I will be the sum of the solutions of two partial 
problems 
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Problem la Problem T in the case 

/(0,r)=*0, 0>7r/2, .. (38) 

but for an arbitrary 
Problem 16 Pioblem I m the case 

. .,(39) 

In lack of strict ladiative eqmhbiium, we have an 

analogous problom Wo confine ourselves to the case (38) 

Problem II I (r, r) ^ 0 and (t) 2 0 aie to be deteiniined 
from (33) and (34), tho incident ladiation being zoio, the net fl.ux 
at tho surface being given, and e/a being given as a function 
of T. 


The Schusle} Schwa) zschild model Tins is the case of mono- 
chroiuatio radiative ocimhbiiiim, with € = 0 The affix v can be 
omitted, the fiocjuency being always the same. On introducing 
an analogous oiitical depth 



tho fundamental ecpiaiions (21), (21') take, according to (26), 

tho form dH-r i) 

„osfl~gtL! = /(T.r)~J'(T,r), ... ,(f0) 


.. .(41) 

According to (22), tho monochromatic not flux is constant. It is 
to be noted that JSddington’s relation (32) holds too For multi- 
Xilioation of (40) by cosd and mtogiation with respect to r yields 

~ J 

If, hero, tho integral is evaluated by means of (41), the integrand 
of the thus obtained ^''-integral is seen to contain tho factor 

Jcos 0y (t, r, r*) dcu = 0. 

Tho vanishing is a consequonoo of (1P0> i o. of the fact that y is 
an even function on the r-siihere, cos 0 being an odd function. 
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Peoblem III Case of infinite optical depth, I (t, r) ^ 0 and. 
J(r,?)g0 are to be determined fiom (40) and (41), the moidont 
radiation being zero and the law of scattering y as well as the flux 
constant F being given 

Peoblem IV Case of fimte optical depth, 0 < t < t* 1 and J' 
are, for a given law of scattering, to bo determined, when tlio 
radiation incident at the ouiei face is zero and when the radiation 
incident at the inner face, 

I (t*, 0 )i 6 < 7r/2, 

IS given. 

In the case of umform scattering, y s 1, Problem III is formally 
the same as Problem la The case y = l and I {r*,0) = oomL, 
9 <TTj2 (black body as a background), is the classical problem 
treated by Schwarzsohild in a fundamental memoir. PiobloniH 
III and IV have hitherto not been discussed in full generality. 
The existence and uniqueness of the solution will bo disoussod m 
Chapter iv Purtheimore, hmits foi the solution will bo givon 
that are independent of the law of scattering T^hoy are the saino 
hmits which were found by Schwaizsohild in the case of uniform 
scatteung In this classical case, considerably narrower limita 
will moreover be found 

§6. FIRST CONSEQUENCES. THE RADIATION 
FROM GREAT DEPTH 

We first diaw some simple oonolusions from the positivity of tho 
mtensity and of the Ergiobigkeit, confining ourselves to the two 
standard oases formulated rn the precedrng seotron. Wo suppoao 

e ^ 0 for all sufflcrently largo t, 

(36) shows then that F (t) never increases, for r suffioionlly largo* 
Prom (36), we have 

K{T)-K{(i)=jyn.)dt. .(.J2) 

Since, according to Z ^ 0, K is never negative, (42) ini]ilies 


( 43 ) 




FUNDAMENTAL PRINCIPLES AND PROBLEMS 19 


This proves that ^ 0 holds for all laige enough r, i.e that the 
not flux IS outward for these r, for the opiiosite assumption would, 
according to the decreasing of Fy imply that the left-hand side 
of (43) tends to — oo as r^oo Those considerations show also that 


the hmit 


F^=limF{T)^0 


exists From this fact and from (42) we obtain the important 
relation 

(44) 


1 C 

Fco = Inn — cos^ 01 (t, i) fZcu. 

T=eo 


The generality of this equation must be emphasized In strict 
radiative equilibrium, 6=0, F<x, is of course the flux constant F, 
(44) holds in the Sohustei-Sohwarzsohild model too, because 
Eddington’s relation was found to hold in that case. 

From (35) wo now get 

7rF^ = 7rFco + 4;7r f ~dr, Fco^O. .(46) 


The second term, taken for t = 0, 

pco g I’M 

4.7T \ - d!r = 47r pedx, 

jo J-M 

lopresents, as it should do, the energy liberated within a normal 
cylinder through the whole slab, having the unit area as its 
cross-section, As a special case of (46), wo note that 

Fco = Fo-4 .(460 

Jo « 

rt is tlius seen that Problem It additivoly decomposes into two 
partial problems, the first one being Problem II in the case 
Feo = 0, while the second one is simply Problem la, F = Fm> 
We theroforo confine our attention to 


Problem II a. Problem II in the case 

i.e in the case where the not flux at the surface equals the energy 
liberated per unit time within a oyhnder of oross-seotion one, 
The solution of the general Problem II is obtained in adding 
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to tho solution of Problem IIa the solution of Problem la with 
an arbitrary F=^Fm'^0 

On solving the equation of transfer (tO) with leaiieot to tho 
intensity, we get 

pr' 

= —sec0 c-^ncoO J 

.. (10) 

It may seem surprising that no boundary condition for t-^co 
appeals in the formulation of Piobleina T, II and III Tho sup¬ 
position, however, that I and J be non-negative quantities, 
allows us to dispense with it 
Theorem I. We have 


/(t, r) = sec^J c )cZi, 0^d<7r/2, 


(47) 


1 e the ladiation coming fiom the deeper layers {0 < 7r/2) iw 
solely due to the Ergiebigkeit of the material within tlie slab. 

Proof Considering in (46) the case 0 <')t/2, t<t', we infor, 
accoiding to J>0, that, for any fixed direction r, ^ < 7r/2, 
e-’-'«°°®/(T',r), O^0<7r/2, 

deoieases with increasing t. According to 7^0 this imphos tho 
existence of the limit 

^(r)=!llm I {r\r) ^ 0, 


T'scitt) 


for 0<7r/2 We are thus allowed to pioceed to tho limit t'~>oo 
in (46), 

7 (t, 1 ) == ^ i {r) -H sec ^ J c- ^ J (i, r) dco, 

for 6 < irfi. The first (non-negative) term on the right-hand sido 
represents the radiation not due to the Ergiebigkeit of tho slab. 
Such radiation cannot, however, exist in the steady state. Wo 
have from (48), according to sec ^ 1, 

7(T,0^e^^(?’)^0, 

thus yielding 

j 7(T,r)oos®07aj^ j 7oos^0dw^e’‘ j iooB^Odo), 

J J+ J + 
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This IS, according to the fimteness of Fco, compatible with (44) 
only if ^ 

IGOS^ddoi — Oi 

V 

10 if ^ (?) vanishes for almost all diieotions ? (in the sense of the 
Lebesgue measure) Values different from zero in a set of measure 
zeio are, however, physically moanmgloss It should be noted 
that wo cannot measure intensities, but only amounts of eneigy. 
Such an amount is always lepiesented by an integral JIdco ex¬ 
tended through a finite cone of dnections The values of J in a sot 
of diioctions of measure zero are heie, of couise, without influence 


§ 7 TITK INTEGRAL EQUATIONS OF THE PROBLEMS 


The reduction to integral equations of boundaiy value problems 
of the theory of radiation goes back to Hilbert and Sohwarzschild 
Milne’s integral aquations of Problems I, II On sotting t = 0 
in (46) and writing t instead of r' afterwards wo find, for 0 > 7r/2, 
the radiation coming from the upper layeis, 

I (t, r) = I1 1 (0, r) + | sec 01 J""e- IsooOfj 

..,,(48) 

When (48) and (47), J~J(t), are inserted into (31), a linear 
integral equation for J (t) is obtained In order to find its kernel 
we intioduoe, according to Sohwarzschild, the functions 






= i j* e-j'|soofl| I aec^"’''^ ] doj. 


. m 


On setting 5 = | seo01 they appear, according to doj = smOdOd{l}f 
in the well-known form 


P CC 


e~^^8~^^d8, 


.(60) 


Insertion of (47) and (48) into (34) yields two essential terms, the 
fh’st one being i ^ 
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due to the energy liberated and to the incident radiation, wlnlo 
the second teiin is the sum of two double integrals 


•Itt 



J {t)dtclo} 





J {t) dtdo), 

.... (51) 


The integrand being measurable and non-negative, we may inter¬ 
change the Older of integration, obtaining thus, acoording to 
(49), 7i== 1, 1 

J{t)J!J{\T~t\)dt, . . (52) 

for (61) The use of a brief symbol A for the linear integral operator 
IS convenient for the treatment of the problem, Wo thus obtain 
IMihie’s first integral eq^uation for the determination of tho 
Ergiebigkeit / (t) = ^ (t), 

J'(t)=A(J), + ^+AI )■)<!«. ...(63) 

A formal disadvantage of (63) is that it does not oontain tho not 
fiux 7 fF{T) An equation containing F is, acoording to Milno, 
obtained by nisei ting (47) and (48) directly into (30) According 
to (49), w = 2, we thus get Milne’s second integral equation 

^(T)=iJ co8ee-^l®''°^IJ(0,r)«^w 

-^2 1* J{t) F2{t'-T)dt'-2 j J (i) F2^{T — t)dtt 

.(54) 

containing all the data given in Problems I, II (54) goes, of 
course, together with (46) 

If (47), (48) are inserted into the integral (32) representing A, 
we obtain, according to (49), n — S, the equation 

A:=i| cos29e-^l»“'’l/(0,r)dM + 2r.;(«)i!8(|T-«|).tt. 

" .. .( 66 ) 
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Since the first term on the iight-hand side is not greater than 

J?’o— QOh^dl{0,))d(jo, 
it follows that (44) can also bo wiitton in the form 

F^^hm~rj{t)E^{\T-l\)dt (60) 

T=co T Jo 

Forthoaotual solution of PiobleinsT, 11,it is moie convenient to 
start fiom the integral equation (63) together with (66) and (46'). 

Theorem II Any non-negativo solution -B(r) of (63), that 
satisfies (60) where is determined by (45'), yields together 
with (47), (48) a solution of Problem 11, with given eneigy 
liboiated, given incident radiation and a given sin face flux 
ttFq > 4:7T^p€dx (integrated through the slab) 

Corollary 1. In case of Pioblein la wo must, in (63), set € = 0 
and / (0, ?) = 0, (? > 7r/2 

Corollary 2. In case of Problem 16 wo must set e = 0 in (63) 
and F—Fat — Om (60). 

Corollary 3 In case of Problem TI a wo must set / (0,0 — 

6 > 'trl^y and = 0 in (66). 

The prool of the theorem is plain from the above. 

The, %nte(j)td equation for P)o()lem III, Since the incident 
radiation is supposed to be zero wo have from (40) 

/(t,)') = [sooflI /((,»■)*: 0>nl!i 

“ ... (67) 

for the radiation coming from the upjier layers. Prom (11'), (11") 
and from (41) we infer that the Ergiebigkoit is always an even 
function of direction, J (r, -r) = </(T,/). Let us introduce the 
lineal integral operator 

A (J) i I J"n (t. »•,(.)') ((, /) dldio', 

. m 

where 11 is defined by 

H (T,r;f,r')s=| soo(l' | (T;r',0.(59) 
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for any two dirootions ?, ? ^ ///Stt is thus tlio koruol of tlio oporatoi* 
A In tho simplost caso of uniform scattoring, )/=>!, this opoiatoi 
leducss to tho simplor Milno-oporator, J boing indopouclont of 
direction For an arbitiaiy law of soattoiiiig, howovoi, no 
simplification is possible Tho meanmg of tho goiioral oporatoi is, 
that, on inserting (47) and (S7) into (41), tho linear integral 

equation J(T,/)=A(/)r., .W 

IS obtained for the determination of tho Ergiobigkoit In analogy 
to Theorem II we have in the piesont case 

Theorem III Any non-nogative solution of (61) satis Tying 
(44), = yields togother with (47) and (67) a solution of 

Problem III. 

The integral equation for Pioblem IV. Since tho inoidoiit ladia- 
tion vanishes, (67) gives the intensity for 0>7tI2 Tho iiitonsity 
of tho radiation from the deeper layers is, howovor, given by 

7 (t, 0 - I (t*, ?) + sec 0 J (/, r) cU, 


0<7tI2 .( 01 ) 

We define, in this case, tho integral operator 

L{j)r,r . m 


by replacing, in the integral in (68), tho iqipor limit oo of integra¬ 
tion by r*> On inserting (67) and (61) into (41), wo obtain tlio 
linear integral equation 

J{t,i)=L (^ f (t*. >')do)' 

.( 83 ) 

for the determination of the Ergiobigkoit, 0 ^ t ^ r*. 

SchwarzschiWs integral equation for umfo) m scattering Smoo, 
in tins case, J is independent of direction, tho operator 
simplifies to i 

J{t)P^{\r--l\)dt, .(04) 

and the integral equation (63) becomes 

J’ (t)=£ (J),+AI «»O' / (t*, e')dm'. 

.( 04 ') 
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Tlie olassioal case treated by Sobwarzsohild is 

J (t*, r) = I* = const, 0<TTj2. .(66) 

In this case the second term in (64') becomes, according to (49), 
= 2, simply ^ 1*^2 (t*-t) . (06) 

For later jmrposes we add the lemark that in the case 

/(t*, 7) = cos0, 6<7Tf2 (67) 

the second term in (64') becomes, according to (49), ?i = 3, 

P3(t*-t) .. .(68) 

Positivity of the lei nel An obvious as well as important pro¬ 
perty of all the integral operators mtioduced above is that their 
kernels take only positive values We thorofoie call A, Tj 2 ^ositive 
operator3. This positivity oooius in all boundary value iiroblems 
of the theory of radiation An evident consequence of that 
property is the 


Lemma For any function (h (i, j') ^ 0 we have A (^I>)T.,r > 0 
eveiywhere, unless 0 vanishes identically. The vanishing is, 
strictly, to be understoo d in the sense of the theory of the Lebesgue 
measure, i o, vanishing up to an inessential (t, ?')-S6t of measure 
zero. 

We may express the lemma also in the following form <I> ^ 0 
and 0^0 implies A(fI>)T,, >0 for all r, 1 . Another equivalent 
formulation is that ^ ^ 


implies the strict inequality 

A (T),,,> A ((!)),., .(09) 

for all r, ir. The equivalence follows fiom the linearity of the 
operator, m particular from 

A(T) = A(a))-i-A(T-(h). 

The same thing is, of course, true for the operator L, 
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§8 SOMB PROPERTIES OF THE FUNCTIONS ^„(^) 

For the discussion of the simpler integial equations we need some 
simple properties of the functions ?i^0 We have 

while becomes logarithmically infinite as 

^^0 Foi 1 all functions (a) aie, however, continuous a-h 
By partial integration in (60), the well-known leouisioii 
foimula = .(70) 

IS obtained. We have, fuithermoie, 

J X 

Fiom (50) the inequality 

*-( 72 ) 

18 obtained. Let us compare the right-hand side of (70) with tho 
right-hand side of the equation obtained from (70) by replacing 
n with w -1 We then find, according to (72), 

.( 73 ) 

for > 1, the equality sign holding only at a’ = 0 We note that 


If the left-hand side of (70) is combined once with (72), anothoi* 
time with (73), we find 


p~X p-X 


. . ..(7B) 


x+n x+n- 

for w S 1 On applying Schwartz’s inequality to 

rw/ fi n-l\ { B n+l \ 

•®n(^) = J (® ^ ® s 2 \(l8^ 

another inequahty 

('^) ^n~\ (^) 1 ('®) ., . , (76) 

IS obtained From (71) and (76) we get 



... . ( 77 ) 
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For later piuposcs the formula 
d E^{x — a) 


< 0 , ' r >«>0 


.(78) 


dx E,, (a.) 

will be useful. For the proof wo note that the left-hand side 
equals, according to (71), 


E^X) 




By means of (77), the quantity within the iiaionthesis is seen to 
be negative 

Some simple integial formulae will, furthermore, be needed 
From (71) and (74), we have 


Partial mtcgiation yields, according to (71), 

1 


4 ^ 




.(70) 

■ -m 


On combining (79) and (80) wo got 

iM, (t - i) * = (t - {) 74 (J) di = M„ n (r) - 

.(81) 


and 


’T* 


'T^f—T 

JQ 




.( 82 ) 





CHAPTER II 

SOLUTION OP PROBLEMS I AND [T 


§9 THE SOLUTION OP PROBLEM la 


Millie’s integral equation of the Problem la is, aocorclmg to (38) 
and (63), e = 0 and J — 

R(t) = A{R)^ .(83) 

Physical reasoning would lead us to the oonjectuie that, in strict 
radiative equihbrium, the radiation becomes nearly isotropic in 
great depth, i e according to (28), that 1 (r ,) )jB (r) ->■ 1 as t->oo 
(44) would then show that B (t) is asymptotically linear, 

B — - fcos^ Odoi . hm , 

TT J T=oo 'T 

i.e. that B (t) = |Pt for laige r. We have, indeed, the 
Thboberi IV. Problem la has a solution 


J {r) = B{r) = \Ff{r)^ .(84) 

where/ satisfies Milne’s integral equation (83) and the inequalities 

/(T) = T + e(T), ^<|?(t)<1 .(86) 


Proof Prom the definition (62) of the operator A and from 
(79),.* = co,,vehav6 i = + . ..(go) 

Purthermoie, from (81) and (82), 1, t* = oo, 


t = A(0^-P3(t). .(87) 

On combining (86) and (87) we got 

T + c = A (^ + c)^ + ^ {c ^2 (t) ~ ^3 (t)} .(88) 

The smallest value of c that makes the second term ^ 0 is, ac¬ 
cording to 9'nd Eg jS ?3 for large r, o = l We thus get, 

for all T, ^ x / A, . 

f{T)>A{f)^, /(t) = t+1. . .. .(89) 

On the other hand, o = ^ is, according to 213^ > , tlie greatest 

value of 0 making the second term ^ 0 for all r. We thus have 

A(t)<A(A)„ f^{r)^T + l .(90) 


for all T > 0. 
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We now set suocessively 

/„+iW=A(/„), ....(Ul) 

The inequality (90) has then the simple form 

/i(t)<A(t), t>0 ...(92) 

Aooording to the positivity of the operator A, we have from (92) 

^ifl)r<Mf,)r> T>0, 

or, with regard to (91), (t) <f^ (t), r ^ 0 Successive application 
of the oimation A generally shows that 

L{r)<fn^i{r), T^O, ...(93) 

1 e each function of the sequence/,^ lies above the picoeding one 
From /i(t)</(t), t^O, 

wo now have, again according to the x^ositivity of A, 
A(r)-A(A).<A(/)., T^O 

From (89) wo got therefore /a(T)</(T). Tf the operation A is 
applied a second time, we get 

fa (t) A (/g)^ < A (/)t , T S 0, 

whence, according to (89),/3 (t)</(t) follows. In tho same way 
wo obtain the general moquahty 

A(t)<A(/),</(t), r^O. .. .(9d=) 

(93) and (94) evidently show that tho limit function 

/(T) = HmA(r) 

Jlra 

exists and that it lies between tho limits A=t + | and/=T-M, 
T^O. Jt follows from well-known theorems that, in (91), we oan 
proceed to tho limit n=‘CO under tho integral sig n. / (t) is therefore 
a solution of/—A(/) with the property (86). The smoothing 
property of an integral operator shows, Anally, that / (t) is a 
oontimious function, 

On sotting S (t) = cf (r), 

it remains to find tho relation between o and tho flux constant F 
From /ssT+iJf and from (81), (82), 3, t*s=co, wo find 

/'« pco 
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The last term being, aooorclmg to 0<g'<l and to (79), 
t'*’=soo, less than f, we infer from (66), J ^cf, J’oo = -?’> that IP 
equals i e that (84) holds, q.o d. 

Prom f—A{f) and fiom (87) we infer that the remainder 
function q (r) satisfies the mtegial equation 

q{r)^A{q),^^]I)M .(96) 

Let us, now, insert (84), f=^r-\-q, mto the second integral 
equation (54), where e and the incident radiation are to be omitted. 
Making use of (81), (82), ?a = 2, t*=oo, wo find that </(t) also 
satisfies the integral equation 

f q{l)Ez{t-T)dt~ f . .(96) 

Jt * Jo 

On inserting (84) finally mto (66), we obtain by moans of similar 
computations, and according to (32), the equation 

.. .(97) 


a being a constant, the determination of which will be post¬ 
poned. 

Piom (22), (34) and (87) we find the temperatuie T as funotion 
of the optical depth t, 

(5f)" = ^’ = UW. .(98) 

Pormulae (47), (67) for the intensity can bo written m the form 

^co 

/(t,?)= e“*J‘(T+soos0,r)iZs; 0<7r/2, .(99) 

Jo 

f rlfiooOl 

e“«J(T+^cos^,r)d5, 61>rr/2. ...(99') 

0 

In the present ease we have / = /(r, 6) and J^J{t)=>B (r) On 
inserting (84), (86) mto (99) we get 




0 ^ 0 < 7r/2, 


.( 100 ) 
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for the radiation coming fioin the deepei layeis, in particular, 
7(0,j^cos 0+ e“®g^(5cos0)d!sj .... (101) 


for the radiation emeiging at the surface 

It IS to bo emphasized that the emergent radiation is in¬ 
dependent of the paitioular behaviour of a within the slab, 
piovided that the material is gray Since Pioblem III for uniform 
scattering foimally coincides with Problem lu, we may also say 
that the emergent radiation (for a given frequency) does not 
depend on the course of ay{x) It should be repeated that (101) 
also holds independent of any assumption about temperature 
In the case of the sun, (101) lepiesents the law of daikenmg on 
the solar disk, in good agrcoinent with the observations, 


§10 UNIQQENESS OE THE SOLUTION 
The nniquonoss of the solution of J^roblem I a follows from the 
Theorem V A non-nogative solution of the integral equation 
{83) IS necessaiily of the foim B (t) == const /(t) 


Pi oof Lot B (t) be a non-nogativo solution of (83) We set 


b = greatest lower bound of 


Sir) 

fir) 


.( 102 ) 


and B*{T)^BiT)~bf{r), 

/(t) being the solution found in the preceding section. B* is then 
also a solution of (83), with the pioporties 

iJ*(T)&0, .(103) 

Or.lb.of5^^ = 0. .. . .(104) 

fir) 


All wo have to prove is that B* (t) identically vanishes. 

Prom (104) and from the definition of the greatest lower bound 
wo infer the existence of an infinite sequence of numbers r„, 
yss 1, 2, 3, such that 


lim 

v=w 


BHr,) 

fir.) 


= 0 . 


(106) 
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Here, we iniist distinguish between two difteient oases, (<?) the 
Ty. have a finite limit point r*, ( 6 ) we have T(,~>oo as v->-c» 

Let ns fiist treat case (a) It is always possible to piok out of 
the Tp a subseq[uence oonveiging to t* Ho generality is lost in 
assuming that as v-^co From (100) and fiom/(T*)>0 

limJS*(r,)-0, 

p=<o 

whence, according to the integial equation, 

limA(-B*)^ =0 

V=scn 

foUmvs. As a function of t, A has tlio property of loAvor 
semioontinuity, the proof of which fact may bo postponed, 
^A(5%^A(j5*),*, 

VP* 00 

This shows with respect to (103) and according to the positivity 
of A, that B* (t) = 0 The lower semioontinuity is pioved in the 
following way. According to (103), we have, for 

T* — S < T < T* + S, 

A (B), 6 i £*'5 (i) m, (t - 1) it + 5 !” S {t) {t - t) it, 

.... (106) 

8 being a positive number less than r*. (In ease t* = 0 tlie first 
term is to be omitted ) The left-hand side is, now, a continuous 
function within (t* — 8 , t* + S), the logarithmic singularity of 
being excluded. We have theiefore 

Jo ^ J 

. .. .(107) 

Since this luequahty holds with an arbitrary 8 , we find for 8 -> 0, 
hmA(.B)^^A(.g)^», 

r~^T* 

1 e. the lower semicontinuity, 

Case ( 0 ) Omitting the star in B* we have, acoording to 
/ (t) r, to prove that a solution B (t) of (83) satisfying 

BiO, hm^^Lo, t.h-co, .(108) 

l/a#c0 '^V 
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vanishes identically Accoiding to and to we 

have c a> 

I^ BE^{\T~t\)(U^2A =2JS(r) ..(109) 

B{r), being a non-negativo solution of (83), determines a 
solution of Problem Xa The ooi’iesponding not flux B is not 
known, however, we know from §6, that nccessaiily F'^Q, 
Piom ( 66 ), J = By and fiom (108), (109) we infer that 

J?’g4lim^^ = 0 , 

V=3CO Ty 

thus yielding F~0 On inseiting this value into (64), J=By 
inoidont iadiation = zoio, we find for t~ 0 

this being compatible with 15 ^ 0 only if B (t) s 0, which com¬ 
pletes tlio iiroof of the uniqueness 

CoROLLABY. Every solution of (83) with a finite lower bound 
has the foim const /(t) 

Proof In the solution B (t) ~|- of (r) the constant c can bo chosen 
so large that this solution is everywhere positive Theorem V 
then completes the proof. 

§ 11. A DIFFER33NTIATION FORMULA 

For the continuation of the discussion of Piobloins I and II, a 
study of the more general operatoi 

K(S)r=f‘j(l)U{\r~l\)dl .( 110 ) 

and of the related integral equation will be useful. Wo introduce 

the notation /»«> 

((/./i)=J^ g{t)h[i)(lL 

The kernel of ( 110 ) is evidently symmetrical The symmetry oan 
also be expressed by the general relation 

to. A (A)) = (A. A (17)) r f “ /i (IT - (I) y (() A (t) Mih. 

.(Ill) 


11 M P 


3 
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Wo start with an important jiroporty of the imrtioular operator 
(110). For this purjioso let us suppose, that II (^) is continuous 
for.-« > 0 and satisfies the conditions 

//(a,) = C»Mog^), x-i-a, U(v) = 0((r% »-i-oo 

.. ..( 112 ) 

Lot h (.t) 1)0 continuous for a ^ 0 and continuously dilferentiablo 
Cor .V > 0. h (a) is, furthormore, supposed to have tho properties 

h ' (r) = 0 ^log 0, h ' { x ) = 0 (e^*“), a oo, 

.... (113) 

with a suitable < 1, Under thoso conditions, tho formula 

™ A {h)^ - A {li\ H- A (0) 1/ (t) . ... (114) 

liolds. 

Pi oof Wo have 

A{A)r- II[i)h{r-t)dt-¥ rII{i)h{T-\-i)dt . (116) 

jo Jo 

On writing (t) for the first integral, we find 

giir gi iiL 11’'' V/ (f) h{r+s-i)dt 

+ .( 116 ') 

Jo o 

Ilore, tho first term obviously tends to A (0) II (t) as 8 0, In tho 

second torm, however, limit process and integration may ho 
intorohangod. Tho proof of this follows the same lino as the pioof 
that a Newtonian potential can bo difierentiatod under tho 
integral sign, oven when the moving point lies within the mass 
That tho second integral in (116) may be differentiated under tho 
integral sign, follows from a well-known theorem, according to 
which this IS always permitted within a certain T-intorval if tho 
T-derivative of tho integrand lies, in absolute value, below a 
fixed intograblo function of t This test apphos hero since, accord¬ 
ing to (112) and (113), 

I II {t) W (t+01 < const, 
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holds for all sufficiently large the constant being independent 
of T. Within an arbitraiy finite r-intervnl, having t = 0 as an 
outside point, the suppositions of the above theoiom are fulfilled 
because the absolute value in c[uestion is less than const 
for large t and less than const. | log t \ foi small t 


§ 12 


ASYMPTOTICALLY LINEAR SOLUTION'S 
OF /=A(/) 

Since the formula of § 11 will not bo used until § 17 we abandon the 
hypothesis (112) about the kernel, Let us consider the integrals 


'00 
J X 


We now suppose that the kernel is positive, 

//(^')> 0 , , ( 110 ) 

and that II{x)^(yIh{x) ...,(117) 

holds with suitable a>0 All those assumptions are fulfilled in 
the case (62), with — (117) implies the existence of //g, 

and automatically of all 11 

The mtogial aquation /(t) = A(/),. .(118) 

Will be studied hereafter by means of l?ourier intogials It is, 
liowever, worth while to consider li first from the sam 0 ,physioally 
more natural, point of view as in § 10. (118) can have solutions 
of various asymptotic behaviour for r large. Under the paitioular 

hypothaaiB /la(0) = i .(110) 

for instance, wo have the 


Theorem VI Under the suppositions made above about the 
kernel, m partioular (119), (118) possesses a solution 

/(T) = T + gf(T), 

whore a is taken from (117) 

Pmo/. In Milne’s case, ll==\Bn (119) is, according to 
// 2 ( 0 )=a j2l3(0)/2=5 certainly fulflllod. Wo first note that ^ 

f (0 ('^) "h h2 ('®).(f 
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Talcing account of (119), we find, after simple oaloulationa, the 
relations 1 = A(1)^ +/^(t) .(121) 

and T = A{t\ — H^{T). .. .(122) 

Accoidiiig to (117), WG obtain fiom (121) and (122) the in*- 
equalities A(t)<A(A).. A = r, 


and 


/(t)^A(/)„ + 


This shows, in the same way as in § 10, that (118) has a solution 
between t and t + l/ty, q e d 
From (118) and (122) wo find the integral equation 

q{r)^A{q), + JJ,{r) .. ..(123) 

for the lemamder function g'(r) We must now remember that 
the above solution/ (t) of (118) was found to bo the limit function 
of a sequence/„ (t) defined by the lecurront relations 

/;n-i('r)===A(/J^, /i = T 

On setting = t + , we find 

dn+i (t) == A {q^)^ + Jig (t) , = 0, 

00 

thus yielding (t) = S A^ (//g)^, 

0 

A”- being the wth iterate of the operator A, A° (/) =/. Tn other 
woids, q{T) IS represented by the Neumann series of the integral 
equation (123), 


Definition If the solution of a linear integral equation is 
represented by the Neumann series, we call it the N-solution of 
that equation. 

The fimotion q (t) of Theoiem VI is accordingly tho N-solution 
of(123) 


Lemma 1. A nowhere negative solution of the inhomogeneous 
equation 


^ (t) = A (^)^ + (t), 0 ^ 0 


.(124) 


18 always the sum of the N-solution of (124) and a nowhere nega¬ 
tive solution of tho homogeneous equation 
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Proof "From (124) we have, aoooi cling to the positivity of A, 

9^ (r) = S A*' (O)^ + A'^+i ((^)^ ^ S A*' (0)^. 

0 0 

Piocecdiiig to the limit n = oo, we see that the existence of a 
non-nogative solution of (124) implies the oouvergenoo of the 
A-seiies of (124), this series consisting of nowhere negative terms 
The lemma follows then from the meq^uality 

f^(T)^2A'^(tI>)^ 

0 


Lemma 2 A continuous solution f/ of (118) with the property 
glf~>c, T-^co, necessarily equals c/ 

Pi'oof, Fiom tho supposition made in the lemma wo readily 
infer that the function ^ 


attains either its maximum or its minimum value foi a fimto 
value of T, say t* Without loss of generality wo can oonflno 
ourselves to the case of the maximum value M, In this case we 

We MJ{r)-a{r)iO, Mf{T*)-a[T*) = 0. 

The left-hand side is hero a solution of (118). According to tho 
positivity of A, however, a non-nogative solution of (118) can 
only vanish somewhoio if it vanishes identically. It follows at 
tho same time that Jlif = c, q o d. 

Lot us now apply those simple results to the equation (121), 
having 1 as a solution All requirements of Lemma 1 are fulflllod 
in this case 1 is therefore nowhere smaller than the A-solution 
of (121) Tho A-soUition being continuous everywlioro, tho 
dihorenoo is an evidently bounded and continuous solution of tho 
homogeneous equation (118). If this solution is identified with the 
g (t) of Lemma 2, / (t) being tho solution of Theorem VI, we find 
from that lemma that g{r) vanishes idontioally. 1 is therefore 
the hT-solution of tho equation (121) Wo collect these facts in tho 


Lemma 3. g{r) of Theorem VI is tho W-solution of (123) 1 is 
the A-solution of (121) 
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§ 13. AN AUXILIARY THEOREM WITH 
APPLICATION TO PROBLEM 11 

Let A, for a moment, Tbe an arbitrary symmetrical operator For 
the solutions of the simultaneous integral equations 

.(126) 

ive find puiely formally 

{<l>i) ^ 2 ) ~ (^1 > ^ 2 ) ~ (•'Ai > A (^ 2 )), (^2 »~ (^3 > ~ (^As > A 

and according to the symmetry of A, 

(<Ai/^2) = (<A2/I>i) .(126) 

This formula is, of course, always oonect when the integration 
interval is finite and when the kernel is oontmuous. The pioof 
fails, however, in all of our coses, since here (<Ai><A 2 ) 1 ® always 
infinite. Nevertheless, formula (126) holds also for certain 
smgular integral equations, of course under certain restiiotions. 

AuxiiirABY Theorem Let A ho a positive and symmetrical 
integral oiieratbr Furthei more, let us suppose that ^ 0, 
0 in (126) Formula (126) then holds, if in (126) <l>i and ^2 he 
the JW-solutions 

Proof, The iterated operators A** aie well known to be sym- 
metiioal once A has this property From 

CO CO 

<^i==SA*'(chi), 

0 0 

we have 

(<ta,^,)=S ($2, A- ((DJ)=S A" (<»2)) = (<I>i, (fa). 

0 0 

The termwise integration of an mflmte senes of functions, applied 
here, is readily seen to be justified, all terms being nowhere 
negative. 

We now return to the special operator (110) and irrove the 

Theorem VII. The necessary and suilfioient oondition that, 
under the suppositions (116), (117), (119), 

^ (t) ~ A (^)y+O (t), (E> ^ 0 


(127) 
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have a nowhoie negative solution (j), is that 

(1,(I))=J 

be finite 

Corollary. Under the same asaumptiona, <I> however having 
any sign, (127) has always a solution if 

(<)!<** 

bo finite The JV-series converges absolutely 

Proof According to (117) and to the fact that H^{%) is a 
dooieasing function, we have 

IJ (i — t) ^ {1 — t)> a/Zg (0> 

thus yielding, with lespect to ^ ^ 0, 

^(0//(i-T)df^« .(128) 

Jr Jr 

Furthormoie, wo have 11 (t — /) § a/fg ('’f' “ 0 since II^ (aj) lies 

between positive limits within (0, t), 

f ^(i)//(T — ^)di^const f <l){t)II^{t)dt, .(129) 

Jo Jo 

the constant being positive and depending continuously upon t 
A (<^)^ must, now, bo finite for some values of r TIonco, and fiom 
(128), (129), we may infer that (^jZ/g) is finite Accoidmg to 
Lemma 1 this holds a foitiori, if <f> is loplaoed by the JV^-solution 
of (127). On applying the auxiliary theorem to the two equations 
(121) and (127) and on taking account of Lemma 3, we find that 

(l,(I)) = (0,/ia)<co. 

The condition is sufficient. On setting 

(;^,^(T) = SA»’(tI))T, 

0 

we infer from the symmetry of A" and from (121) that 

(<)t„, Z/g) = S (A>^ (O), ZZg) - S ((!>, A-' (Z/g)) 
«(q&,SA''(ZZg))g((I>,l), 
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proving thus that {(f)II 2 ) lies below a fimto limit indeponclent 
of n The fact that incieasing set of funotions 

^ ^n-hl > W — 1, 2, 3, 

implies the existence of the limit function (f) (r) — Inn (t), and 
womfeithat (^.if,) = hm(f,./4) 

18 a finite quantity At the same time wo readily infer that 

(^,/4) = (aM) 

The corollaiy simply follows from the theorem, since an 
absolutely integrable function <5 (t) can be written in the form 
^ 1^1 ^ 1^2 both non-negative and integrable. 

Application to Pioblem II Milne’s integral equation (63) is 111 
the case of Pioblem II, J~J3, 


Bir)^A{B), + ^ .(130) 


Supposing that everywhere e ^ 0, we find from Theorem VTI that 
the necessary and sufficient condition that (130) have a non- 
negative solution is, that 

47r|^ ^-dt, .(131) 

i.e. the eneigy liberated per unit time in a normal column of 
cross-section one, be finite We prove, furthermoie, the 

Theorem VIII Suppose that e ^ 0 everywhere The necessaiy 
and sufficient oondition that Pioblem II« has a solution is that 
(131) be finite The uniquely determined solution is given by the 
iV-solution of (130) 

Pi oof, Only the assertion oonceimng the iV-solution of (130) 
needs to bo proved We remember that Problem Xia is the special 
case of Problem II whore the surface fiux tt J’o satisfies 

. 

It is therefore to be proved that the A'-solution B{r) of (130) 
automatioaliy satisfies (132). Prom (64), t«0, J = J3, taken 
without incident radiation, we find 
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B being the iV-solution of (130), we get from the auxiliary tlieoiem, 
applied to (121) and (330), 

which completes the proof of (132), IToi any othei non-negative 
solution B of (130), howevoi, the left-hand side of (132) would 
be gieatei than the right-hand side, 

^-dt ....(133) 

J 0 

In 01 del to obtain the solution of the general Problem II, wo must 
lomember foimula (45') giving when Fq is given 

Theorem IX Suppose that € ^ 0 everywhere The solution of 
Problem II is uniquely detoimined by 

B{r) = B{r)-VjFJ{r\ 

Fro being determined from (46'), B being the .^-solution of (130), 
and/(T) being the function of Theorem IV 

The uniqueness of the solution follows easily from the fact that 
any non-negative solution of (130) is the sum of the jV-solutioii 
of (130) and a non-negative solution of the homogeneous equa¬ 
tion, the latter being oharaoterized by Theorem V. If vF denotes 
the flux 111 Problem TIa, we have 

F^-F^-l-^ Fro> f ... ..(134) 

If the assumption c^O (everywhere) is abandoned, similar 
theorems hold. Put wo have then to formulate the general con¬ 
ditions undo)’ which the solution becomes non-negative as 
loquirod in the physical problem 

§ 14. OTIfEB APPLICATIONS 
THE BOUNDARY TEMPEBATXJRE 

We set, j,5(T)='L(l±^z£W, .( 135 ) 

where /(t) is the solution of (118) mentioned in Theorem VI, and 

] ^t+8 

i?8(T)«|J^ IJ{t)f{r + h-i)dt, .(130) 
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8 being a positive number From the integral equation (118) and 
from the identity (116), we find 

^/8(T)=A(^7s),+ G^5(r) . . .(137) 

Proceeding to the limit 8 0 we get, puioly formally, the integral 

equation foi the derivative/' (t) However, as we do not know if 
this derivative exists, we had better use (137) as it stands 
In order to piove that is the JV-solution of (137), we fiist 
note that n 

!ir8(T) = SA'’(G8),+A“+M!/s)r -(138) 

0 

Keej)iiig 8 fixed we infer from Theorem VI that (r) is a bounded 
funotion of T, | gs I < whence 

|A"+i(9'8)rl<C^A’Hi(i)^. 

According to Lemma 3, we have 

A'^+i(l)-|S A‘'(J2;,), 

showing thus that the remainder in (138) tends to zero as «.-^co. 
flfg (r) is therefore the -V-solution of (137). 

Theorem X The value of the solution /(t) of (118), given in 
Theorem VI, is determined by a simple formula, 

/,=V257(0) 

PiooJ Applying the auxihary theorem of §13 to (121) and 
(123),™ obtain = = .(13B) 

Since, by partial integration, we find altogether, 

according to /=T+<?, = 2^4(0). .... (139') 

Since, furthermore, g^ is the iV-solution of (137), the auxiliary 
theorem can be apphed to (123) and (137), yielding thus 

((Z,G^8) = (i7a,^4) .(1^0) 

On the other hand, we find by means of (136) 

(SiJh)=-l //Wff. («)*+J”/(«) 

The fraction in the second integral hes between the limits 
Hg {t — S), //g (i) 
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because 11^ — is a decieasing function We may tlieiefore 
proceed to the limit under the integral sign, finding thus 

lim(!/s.//,)-=-/o/4(0) + (/./4) .. (Ul) 

0 = 0 

Furthermore, wo have from (136) 

hmC?s(r)=/o7/(T) 

8=«0 

On inserting the right-hand intogial in (136) into the left-hand 
side of (140) it IS seen without difiiculty that we may proceed to 
the limit S 0 under the integial sign on the left-hand side of 

8«=*0 

Finally, wo have from (123), t = 0 , according to g'o=/o, 

/o = (<?,//)+ 74(0) .(1A3) 

Combination of (140), (141), (142) and (143) loads to 

whence, according to (139'), the required formula folloAvs 

Lemma 4 Lot 77 (^)(^) bo bounded for x laige enough. 
Under the condition (112), the solution /(t) of (118) has then a 
positive derivative /' (r), being the ^-solution of the mtogial 
equation /'(t)=A(/'),+/o//(t) .. .(!«) 

Proof. First we note that this applies to Milne’s case, H = 
since as jr->-oo. Since the lemma is only needed for the 

proof of Theorem XI wo shall content ourselves with brief in¬ 
dications of the rather lengthy proof. The mam part of the proof 
consists in establishing an inequality for the jW'-solution of the 
integral equation 

(7 (t) = A((7),-1-477(t). . .(146) 

The proof is complioated by the possibility of an infinite 77 ( 0 ) 
(this is just the case with Milne’s equation) It is tlioreforo oon- 
vonient to split II (t) into two i>arts 

/qU .(140) 

(477, 0 <T<a, 
lO, T^Cf. 


where 




.(147) 
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03 ( 7 ) 18 zero for small t and coincides with /qH (t) for largo r. 
We consider then the two equations 

^i = A(^i) + Oi, 9!»2 = A (^ 2 ) + • .(148) 

separately Accoidiiig to the supposition made in the lemma we 
have (Da (t) < const (t) fox all r, The iV-solution of the second 
equation exists therefore as seen by comparing it with the 
equation (121) At the same time we find 
0 < ^a (^) < const, 

for all r As to the fiist equation, we begin with considering the 
expiession 


J 0 


r>a 


We note that A(<Ih) is a continuous function of t This holds, 
accoidmg to (112), also at t = 0 Foi t large we find, according to 
the hypothesis made, and since deoieases, 

A (Oi)^ < const. J //a {r~t)E{l)dt < const //g (t — a),. (149) 
the constant being independent of a Now we have, for t large, 


^(a(r-a)~Ha(r) = 


rr 


Udt 


r~a 

< 


const, n^dt < a, const. J/g (t — a) 

J r~a 


Since here the constant is again indejiendent of a, we may choose 
once for all a so small that a const < For this value of a, we 
flndtW /4(r-«)<2/4W, 

thus yielding, with lespect to (149), 

A(<I)i)^<const Ha(T) 

for T large The left-hand side being continuous for all r, we can 
choose here the constant such that the inequality holds for all 
T ^ 0, Comparing again with (121), wo find that the W-series 

SAMOj), 

converges for all t, On adding here the term v— 0, namely Oj, wo 
obtain the complete N-solution (j)^ of the flist equation (148). 
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Altogether we see that the iV-solution (j{r) of (145) satisfies an 
inequality (t)< if (t) + ooiiet 

From (137) and (146) we now find 

(7S"“(/ = A((78 —</) + (G'§ SqH) 

The second teim on the right is easily seen to have the following 
properties as 8 -?■ 0. It is 0 ^log-j fox all small t, nmfoimly with 
respect to S, while 

n,{r) 

tends to zeio uniformly fox all t S tq , Tq being an arbitrary number 
greater than zero, A similar splitting leads then to the result that 
(j is the limit of f/8 ,1 e that g (t) ~J' (t) At the same time wo find 
an inequality /' (t) < ii (t) + const ..(160) 

The above oonsiderations show that/' is the iV-solution of (144). 

Theorem XT Suppose that tI)(T) is continuous at t=0 
Under the hypothesis made in Lemma 4, the iV-solution of the 
integral ectuation (127) has, at t = 0, the value 

.( 101 ) 

Jo 


P)oof Since, under the stiictor hypothesis about II mado m 
Lemma 4, the derivative /' (t) exists, wo can dispense Avith the 
complicated formulae (136), (136), (137) and apply the auxiliary 
theorem of § 13 directly to (127) and (144), getting thus 

From (127), T==0, wo find (<^,^/) ——®o> whence (161) follows, 
The reader will notice that (161) can also bo written m the form 
of a Stieltjos-intogral 

.( 162 ) 

Jo JO 

It is, of Qourse, easy to obtain Theorem X from (162). Applica¬ 
tion to (123) yields, according to (Zo=/o> 

which together with (139') proves the theorem. 
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Application to the boundary tempemtioie in Problems I and II 
Ill th.6 case of l^roblein lit we have 

2//,(0) = ^,(0)-i 

. „ 1 _ 


yielding thus 


. ... (163) 


for the solution mentioned in Theorem IV According to (84) wo 


tlieiefore get 




.(164) 


Under local tlieimodynamical equihbxnim we thus flndj according 
to (98), the value 

m A V “ m <i 


foi the surface temjierature of model la, independent of tho 
variation of the absoiption coefficient within the slab 
Let us, now, find a formula for the boundary teinperatuie of 
the model II a. In order to apply Theorem IX we first suppose 
that 6^0 everywheie We find for the A-solution P (r) of (130) 

A-- r-fdt=j [”-*+)■ r-q'dt. 

^ a /oJo® /oJo<5C fojQ^y 
Aocoiding to (132) and (163), the first term on the right equals 


From Theorem IX we therefore obtain, on account of Fq = Fq + F^, 

.(I5e) 

a 0 Jo'* 

for the determination of the surface temperature of model II. 

It should be noticed that (166) holds also in the general case, 
€$0, provided that 

lo 

integrated through the slab, is finite In the next section we shall 
piove that g'(r) increases with increasing r. Anticipating this 
result, we infer from (166) that, in the case e ^ 0 everywhere, tho 
boundary temperature is greater than in strict radiative equili¬ 
brium, 6=0, the surface flux ttFq being prescribed, 
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This minimal property of the boundary temperature in the 
case of strict radiative equilibimm implies, conversely, the in¬ 
creasing of tho lemainder function g(r). Otheiwiae the onoigy 
libeiated could evidently bo distributed in such a way that that 
property becomes violated 


§16 PUOOF THAT ? (t) IN'CREASES 

Let us, for a moment, leconsider the hypotheses made in the 
general case tieated in §§ 11-14 Tho hypothesis (112) concerning 
the behaviour of 11 (a;) for small a was only required for tlio proof 
of foimula (114), of Lemma 4 and of Theorem XI (existence of 
fir)) 

For tho proof of the Theorems VI-X we needed, however, only 
the suppositions stated in (116), (117) and (119), tho latter being 
of particular importance because it is rosponsiblo for the asymp¬ 
totically linear chaiactor of the solution, 

Let us, now, throughout this section, suppose that II (a*) is of 
the form 

II{x)^\ .( 160 ) 


p (s) increasing with increasing s, (119) is satisfied if and only if 


From 


rini^) 


j: 


' dp (s) 


1 

'2* 


-sj: ol-Ti 


dp (a), 


(1600 


wo see that (117) is certainly fulflllod, with cf = I Tho last integral 
IS readily soon to converge also for n^O, making thus tho formula 
valid for all | 0, 

On applying Schwartz’s inequality to II ,in tho same way as 
we have done in deriving (76), wo find 

II,,^<IIn-Jfni 

for all implying thus 

d 

dx IIM 


( 167 ) 
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The analogue of (78) subsists tlierofore m the present goneial 
case, the proof being the same, 




< 0 , %>a 


... .( 168 ) 


We now begin witli the proof that, under the hypothesis (166), 
(166'), the remamdei function (r) of Theorem VI inoieases with 
incieasing t 

Transfoiming the integral operator A by partial integration, 
we find the identity 


« (t) - A {u\ s UqH^ (t) H- (t ~ 1) du {t) — //a {i - r) d-u (4), 

.(160) 

I Stieltjes-integrals being employed beoause wo do not wish to 
make any use of the differentiability of/(T). (123) can therefore 
be written 


J^// 2 (T-i)dgr(i)-£ IIz(t-r}dq(l)^IIa(T)-^oIl2(T) ...(160) 

Let, in the sequel, u(t) denote any positive and nowheio 
decreasing function that satisfies the integral inequality 

u(r}^A(u}^-hIIs(r) .(161) 

for all T. Continuity of u is not required An example of afunotion 
u (t) is furnished by ^ == 1 , for in ( 121 ) wo have > 11^ The sot 
of values taken by all functions u (t) at a particular point t has 
a greatest lower bound U=U(r) u^O and ( 101 ) implies 
17 (t) ^ 7/3 (r) > 0 Furtheimore, C/(r) is evidently nowhere 
decreasing From (161) and from 1/ wo see that 
u(r)^A(Ul^^TI,(r) 

holds for all funotions u, and Uieiofore for their greatest lower 
bound £7 U (t) is thus itself a funotion u, in fact, the smallest 
function of that class. 


Lbmma U is continuous for t ^ 0 Jn tho inequality 

£7 (t) ^ A (‘U.)^-h 7/3 (t) .( 162 ) 

the equality sign takes place at all points not being inner points 
of a oonstanoy interval of £7 t == 0 is included horoin. 





SOLUTION OF PROBLEMS I AND II 


49 


We posi^jone the lengthy proof of this main lemma to the end of 
the seotion and show first that TJ oannot have a oonstaiioy in¬ 
terval, theieby proving that U moreasea and satisfies the same 
integial equation as q{T) This would, liowevei, imply U 
since a bounded solution of the homogeneous equation vanishes, 
according to Lemma 2 Supposing that U have a constancy 
inteival wo find, according to the lemma, from (1S9) 

and (103) 


TJ 


(r) 






Jo IIM J- 


*//„ it-r) 


dU 


= 0, T = c^, 
, ^ 0, T> tf. 


. (163) 

T lying between a and ^ This can, howevei, be shown to lead to 
contradiction Aocoidmg to (167), the second term decreases 
with inoieasing t According to (168), and to dU^O, the tliiid 
teim certainly never increases. TTuitheimoxo, — 
incioases with inoieasing t, the denominator dooicasing and the 
numerator increasing The foiuth term in (163) is thorefoie also 
tt nowhoio increasing function of r Altogothoi wo find that the 
whole left-hand side of (163) decreases, a<T</S, in obvious 
contradiction to the ilght-hand side. 


P}OoJ of the lemma T'or a nowheio decreasing function TJ the 
limits Z7(t-- 0) and U (t*|- 0) exist always. Since V is the smallest 
of all functions avo find 


V{T^Q)^U{r), 

since otheiwiso u (r) = £7 (r — 0) ^ 17 (t) would be a smaller func¬ 
tion u, 

Tn Older to prove the contmiuty of U wo must show that 
t; = 0, whore 37 ; = £/ (tq -I- 0) - £/ (tq) ^ 0.(164) 

The right-liand side of (102) being a continuous function, we 
derive from (/ (t)'^U (tq -t- 0), t > tq , 

£7(t)^A(££),-!- 77a(r)-I- 7 ?, ro<T<ro + S, .(166) 

with a suitable S > 0. On introducing the auxiliary function 


, , ^ ff, r m (tcToH'S), 

(0, T not in (tojTo-HS), 


.( 100 ) 


II M V 


[ 
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we filst find the obvious inequalities 

A (/i)^ -b 1, T in (to , To + S), 
A(/i)^, T not ill (to, To + 8) 
It 18 geometrically evident that the fuiiotion. 


h (r) < 


.. .( 167 ) 


-ji (t) = Cl (t) — 17 / 1 - (t) .(108) 

nevei decreases and that ^^^0 ]?iirtheimor 6 , the second in- 
©quahty (167) shows that, outside the interval (ro,ro-l-8), u 
satisfies (161). On the other hand, we infer from (166) and from 
the fiist inequaUty (167) that (101) is also satisfied by u within 
that interval u is therefoio a meinbei of the class intioduced 
above, and which obviously implies 17 = 0 

It is immediately seen that, in (162), the equality sign must 
hold at T=0 Otherwise we could make the value of V (0) slightly 
smaller without affecting (162), in contradiction to the definition 


of the smallest w-function 

A point T shall be called a ‘proper point’ of U[r) when 
TJ (tO > CC (t) holds foi every t' > t According to the continuity 
of C 7 , an ‘improper point’ is thus either a left end point or an 
inner point of a constancy interval of CC A point that does not> 
he inside of such an intoi val is therefoie oithor proper or a left end 
point of such an mtei val In the latter case, and for r > 0 , the point 
IS oeitamly a right-hand limit point of proper jioints. According to 
thecontinmtyof XJ it is thus sulfioientto prove the equality sign in 
(162) for all proper points tq . Wo define a number 17 ^ 0 by sotting 

217 - CC (To) - A (CC),, - Ih K) .(109) 

and show that '>7 > 0 is impossible As U is continuous, (166) must 
hold in a sufficiently small mtoival (t^, tq-I- 8 ). This time we iiso 
another auxiliary function 

7 ./ (t)- CC(to), t in (To,To-hS), 

1 0 , T not m (tq , To H- 8 ). 

According to 0 < CC ^ 1 wc have 17 < 1, (t) is again soon to 

obey the inequahties (167). The function w, defined by (108), 
belongs then again to the above olass, since 11 satisfies (161) and. 
decreases nowhere (geometrical evidence) This completes the 
proof of the lemma. 
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§16 OTHER PROPERTIES OF (t) 


The eq[uation (96) oan be immediately generalized, 

jy (t) (t -t) - [Jj it)Ih{T-t)dt=SM ■.(170) 

This equation oan be ohecked by diiect differentiation, taking 
aooount of (123), There should be an additional integration con¬ 
stant on the right. On comparing, foi t = 0, with (139) we see, 
however, that it has the value zero The left-hand side of (170) 
oan easily be transformed by partial integration, thus yielding 


J"fl-3(|T-<|)rf!?(i) = H,(T)-2o//o(T). . . .(170') 


From (167) and from Jl,^ follows the existence of the limit 

A = liini-^^1, 

x=^co-^^n (-I 


thus yielding = where 8 = 8 (,t)-» 0 as aj->oo. 

According to (166) we have 

/»1-{C /’1 + 6 

Hn i-i (^t) > J e~^^ 3 ~^’'dp (5) > J s-^^^dp (5), 

p increasing and e being arbitrary, whence follows A^l. This 
implies A= 1, 


From (168) we infer, furthoimore, the existence of 

a!»oo 

From —//,/ wo obtain, by diCterontiation of this 

limit relation with respect to y, (y) = (j> {y) This implies, together 
with 9^(0) — !, ^(y) = ei/. According to (167) and (168), we find 
at the same time that 

/4(T~i)>e^J4(r), r>t. 

Since, now, q{t) increases, wo find fiom (170') 






i. 


e^dq{i)^l-qQ. 


thus yielding 


.(171) 



62 


SOLUTION OTP PEOBLEMS I ANB 11 


Let US, m (166), make an additional hypothesis concerning the 


behaviour of p (s) at s = 1 We suppose that 


(“'^'>(*> = 00. 

Jl 5-1 

... (172) 

When p has a continuous deiivativo this is certainly true if 
p'(0)>0 This takes place, for instance, in Milne’s case, 
p{s) = l log 8 Under this hypothesis we always have 

j^e^dq{t) = l-qo. 

....(173) 

It 18 remarkable that this general relation does not contain 77 or 
p explicitly. 

P?oo/ of (173), We first show that the functions 


^71 (^)—( 5 ) ^^ds ,. 

... (174) 

have the properties 

2,,(r)~>oo, r-^oo. 

.. (174') 


According to as s->co, the second property is seen 

to be a oonsoc[uonc6 of the first one. Interchanging the order of 
integration we find 

^00 roool—71 

~dp{3), 

which IS infinite according to (172) 

Wo need the inequality (c = 7/4 (0)) 

+ . m") 

0 ® ' ‘ i CCS i>x, ^ ’ 

Por t > X the left-hand side equals, indeed, 


'T! 

( 

JO 


e'^Z/g (i ~-r)dr<e^ j II ^ {I ~r)dr< ce®' < ceK 
Por iga:, howevoi, it equals 

t z' cc 

6^773(^—7)77+ e'^ 77 g (t —i) 7 r. 

JO J< 
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The first term, is smaller than oc“^, while the second teim becomes, 
on substituting T=:f + s, 

11 ^ (s) ds, 

being obviously less than (^), which proves (174") According 
to (174") we get, on multiplying (170') by and integrating 
through 0 < T < ^, 

/»» pco 

k ('«) ^ok ('«) < ft (^) + c} e.^dq-\- G^dq 

Jo J X 

Dividing ihiough by l^{x) and proceeding to the limit oj-^-co, 
leads, according to (174'), to 

^00 

l-q^^^^ddq, 

whence, with respect to (171), (173) follows. 

By paitial integration, (173) can be transformed into 


j. 




.. .(173') 


NumG) ical remaiica on Milne's case, Let us first find the value 
of the constant a m the formula (97) The limit of the left-hand 
side as T~>-oo evidently equals the limit of 


^ < 
J U 


q{t) ]UQ{\r-t\)dt, 


Cl) being an arbitrarily fixed quantity. This integral is, now, 
included between the limits 

q{oi) \ B^{\T-t\)dly f B^{\r-t\)dL 
J OJ J (li 

Here, tho integral tends to § as t->oo, whonoo 

§g(ci))^a^f(7„, 

and according to tho arbitrariness of co, 

u=§(7«, . .. .(176) 

TTrom (97), t = 0, and fiom (176), 

{q.B,)^lq^~h .(176) 

Furthermore, from (96), t —0, or fiom (139), 

.( 177 ) 
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Practical computation of a number means inoluding it between, 
sufficiently nairow limits This can sometimes be done by means 
of mathematical artifices. The two lelations (170), (177), to* 
gother with the fact that q inoreasos, suffice to find the value of 
up to 1 per cent, (178) 

or, in rounding off the second digit, = 0*71. 

We have, in fact, from (176) and (177), 

(Zoo = I ((Z> ^3 - + h 

which, aecoiding to 2.E?g > JS^ for t > 0, implie.s 

(Zco < |(Zoo (1> •^3~ i-^2) "Hi” i(Z«> +1, 
i e the right-hand inequality of (178) On the other hand, fronii 

3^?4>2J?3, ajq 

{}J},-Wdq>0 

Since here the integrand vanishes at Z = 0 and at Z = oo, wo obtaiia 
by partial integration < 0. 

Together with (176) and (177), this implies the second inequality 
of (178) 

We now considei the law of darkening (101) and notice that the 
second term in tho parenthesis dooreasos with incieasmg 0, 
0< 7t/2. Milne’s first approximation consists in replacing (Z by uu 
constant such that the emergent radiation gives the required nob 
flux ttjP, thus yielding the value f for tho constant F being a 
definite average of the omeigent radiation, wo see that tho 
second term lies, for 6 near rrl2, below Milne’s approximation % 
and rises, as 5 deoroases, above it Tho smallest value of the second 
term equals = l/-\/3, while the greatest vahio ((/, e"0 surpasses f 
We may add that 3<(5,H)<O.60. .(J79) 

We have, in fact, from (177) and from 

(2, = ((Z. 6"'- -Ba) +1 < ^ -h i < 0*69. 
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§ 17 PROBLEM T& FOR PARALLEL 
INCIDENT RADIATION 

Parallel radiation oanying finite energy is to be consideied as a 
limiting case. Let 

= d'<l, 

be the diiection of an incident lay, and let Aw bo a solid angle of 
directions containing r' We sot 

(7r>S/Aw, r in Aw, 

0, r not in Aw, 

for the incident lacUation. When Aw shiinks down to the single 
direction r', we speak of paiallel radiation of the diioction r' and 
of the flux 

through the unit aiea normal to r\ On proceeding to the limit 
Aw 0, Milne’s fundamental intogial equation (63) becomes 

i3(T)==A(B)^ + |e-^ s = scoO'. . ..(180) 


/(0,t) = 


We define (180) os the integral ocxuation for paiallel incident 
radiation of direction 

According to Theorem VII, (180) has a positive solution. Let 
ns, in particular, considei the A'-solution, i.e the smallest 
positive solution of (180). Wo put 

^(r) = %(r), ..(181) 

l/g being the V-solution of (180) for S~1 

Theorem XII. The V-solution of (180) can be explicitly ex¬ 
pressed in terms of the solution / (t) ^r + q{r) of Problem I a, 

a . (t) = K (/. «-•') 1 '^’ .(182) 

PiooJ. (182) can also be written m the form 

9 b ('t) =* f (/»1 lo ^ ■ 


(182') 
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First W6 sliow tliQit tho function 


ll{r 


-J. 


Mt-r 


dt 


( 182 ") 


solves (180) with a suitable S Fiom (182"), 

A(0) = 0. .. (183) 

Formula (114) is applicable, h being continuously diffeientiablo 
and satisfying W = 0 (t) for r large On sotting 
p (t) = (t) — A (70tj 
we find from (114) and from (183), 

8p + p* =/“ A (/) = 0 

Hence p (t) = const e-®^. Now we show that in 

gs{T)=A{gX+ae-^^ , (184) 

a equals J It is obvious fiom (182') that g^ (t) is a bounded funo- 
tioii of T The boundedness clearly implies that it is the iV-solution 
of (184) Applying, now, formula (152) to (181), O = ae-«^, wo find 

/o 

If j/g (0) IS computed from (182), we obtain, according to/o = 1/V3, 
a= J. This completes the proof of our theorom 
We note the special coses t = 0, oo, in (182), 

!7.(0) = t«• . ■(180) 


thus yielding 


.(1S6) 

Oa (co) V3 

(185) being expressed in terms of the solution 7'of the Problem la. 
In particular we have the formula 


_socf?' 


.(187) 


for the ratio of the limiting iemporaturos 
The fundamental limit relation (66), J^B, Faa — B, is easily 
seen to remain valid m the case of parallel incident radiation, 
Aooording to the boundedness of J5 (r) tho right-hand integral 
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leprosents also a bounded function of t, whence F = Q We may 
collect the above lesults in the 

Theorem XIII The ^-solution (t) = (t) of (180) corie- 

sponds to vanishing net flux, F = 0,io it gives the solution of the 
Problem 16 (insolation pioblem) for parallel incident radiation of 
direction ?' = (tt — s = sec 6' 

By partial integration we find fiom (182'), acooiding to (185), 
that Qs (t) - ffg (co) equals 

A(T)=re»<<ig(()-(l-!/„) + l-i 
Jo 

up to a positive factor For 1, we have, according to (173), 
A < 0 for all r, wheieas foi 5 > 1, 

A(0)=^3-i. A(oo)>l-i>0. 

A(t) being an increasing function, we infer that the equation 
A(t) —0 has at moat one root Furthormoie, on dilTerentiating 
(182'), and using (173), we find that, for s ^ 1, (t) increases with 

increasing r. Altogether we note the following bohavioiu of the 
solution jB (t) of (180); 

s g 1 • B{r) increases, 

v/3 > 5 > 1. B{t) lies first below the value B^o and then rises 
above it 

5 > VS * B{r) lies ontiioly above B^ ♦ 

5^ 1(5 = 1 corresponds to normal mcidenoe) is the only case 
where B is inonotonio 

Oeneralization. Model 16 with parallel incident radiation, and 
with F=0, mainly applies to the upper layers of a planetary 
atmosphere, The hotter solar radiation will actually have a 
smaller absorption coefficient than the cooler radiation of the 
atmosphere The next idealizing step would thus consist in 
making the former coefficient a constant fraction n oi the latter 
ooeflfioient a. 



58 SOLUTION OF PROBLBJrS I AND H 

In Older to find the integial equation we must go back to tho 
fundamental flux equation (64) The flist teim on the right is tho 
only one that contains tho radiation from outside In this term, 
theiefoie, t must be replaced by nr Proceeding as bofoie to tho 
limiting case of parallel radiation, we find, J ~ B, 

F = Sco3d'e-^^^^<^^(‘'+2j BE 2 {i-T)dt- 2 j^ BE^{r-t)(ll, 

1^=const On differentiating this i elation and on dividing 
thiough by 4, we find 

Since this is of the same type as (ISO), all formulae of this section 
remain unchanged, with the only differenco that S is to bo 
replaced byw^S and that s is now sec 0' We have, instead of (181) 
and (182), 


B 

5 (t) (t)) 


sec O' 


§18. THE EMERGENT LIGHT 
It wiU be convenient now to mtioduce tho now notation 

rg(ff) = ff 1^ 0(cr) = arj ^ .,,.(188) 

The emergent radiation (law of darkening) becomes then 

Problem la* I (0, 6) = ^m (sec d), .(189) 

and in the case of parallel incidence, 

Pioblemlfi* 7(0,d) =/Sireoofl'(seed) ... .(189') 

Prom (183) we find, by partial integration, 


(‘to M 

er<^ni'dt = a\ 

Jo Jo 




thus yielding, again with respect to (183), 
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This gives, together with (182) and (182"), the formula 

4 5 - 1 - 0 ' 

We have theiefore the 


Thfoebim XIV In Problem 16 for paiallel radiation with the 
angle of incidence 6\ the emergent radiation can be exphcitly 
expressed in terms of the same radiation foi Problem la, 


I(O,0) = f^ 


cos (sec 6) 0 (see O') 
cos 0 + cos O' 


(191) 


In the general Problem I, with parallel incident radiation, but 
with an arbiti ary flux ttF ^ 0, the emergent radiation is, of course, 


I{0,0)~- 


F + S 


cos (see 0) 

--~ 0 (sec 6 ). 

cos 0-1-cos 6 


This applies to binary stais, whore the incident radiation is due 
to the other component, 


§ 10 PROBLEM lb FOR ARBITRARY 
INCIDENT RADIATION 

Since arbitrary incident radiation can bo obtained by super¬ 
position of parallel bundles of difleront directions, wo should bo 
able to express the solution of Problem T 6, for any given incident 
radiation, in terms of the solution for iiaxaUel incidence, Por 
given incident radiation J (0, r') > 0, a solution of(63),e=i0, J = 

IS, in fact, obviously given by 

5(T)=i J_/(0.r')cr.(T)(««,' 

TT 

1 r 27r /* 2 

= - -f(O,7r-0',f )0BeoO (T)sin0'd0'#'.(192) 

'rrJo jo 

Let us sliow that this solution corresponds to vanishing net flux, 
JP = 0 We note that, for any value of t, (t) decreases with in¬ 
creasing 3 This IS physically evident, for it means that the 
temperature becomes smaller when the incident radiation is more 
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oblic[ue Remember, in fact, that (t) is the iV-solution of (180), 
s = 1 The above assertion follows simply from the fact that tlio 
second term on the right of (180) deoieases\vith inoieasing s, 5 ^ 1. 
We have, in particular, (jg (t) ^ g ^ (t), and fiom (192) 

The boundedness of B (r) implies, as before, the vanishing of tlie 
net flux. 



CHAPTER III 

DISCUSSION OE PROBLEMS III AND IV 
^20 SOLUTION OF PROBLEJI III 


Pioblein III has been analytically sot up m § 7 We remember 
that the flux ttF is constant and that the limit relation (44) holds, 
in consequence of the liypotheses, that the law of scattering is 
symmetric in the directions of the incident and scattered lay 
(11'), and that the same amount of radiation is scattered in 
opposite directions (11"). 


Tiieoubm XV. Problem III always has a solution 

«/(T,r)==fE/(r,r), .. .(193) 

with the property 

/(T,0 = T + {Z(r,/), 0<(?<1, ..(194) 

the limits thus being independent of the law of scattering. 
Under the hypothesis 

y(T,r,/)> const >0 

the solution is unique 


P)ooJ The reasoning is essentially the same as in Problem la, 
1 6 in the special case of uniform scattering, y ^ 1 Wo put 

J‘j^(T,r) = T, i/(T,r)=T+l, .(196) 

and denote by the rospcotivo intensities found from (47) 

and (67) in sotting / = J = On inserting these intensities 
into the right-hand side of (41) wo obtain the expressions 

respectively. 

Using the more convenient relations (99), (99'), instead of (47) 
and (67), we readily find 


/ / \_„ (?<7r/2, 

^) == cos0-h I cosI 0>7r/2, 

1 7/ ^ fT4-l + cos0, 0<7r/2, 

and U»U = |^^;L + cos 0 >-(l + cos0)e“Tisof'fll, 0>7r/2 

These equations imply, in particular, the important inequalities 
(t, 7 ) ^ T + cos 0, / (t, 1 ) ^ T -j- 1 ■+■ cos 0, 
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for all directions ? — (0, Inserting this into (-ll) and taking into 
account the relation 

J cos d'y (r, ??’) dco' = 0, 
we obtain the important inequalities 

A{J)^^r<T + l^J (r,?). 

According to the positivity of the operator A, we oan, now, apply 
the same leasomng as in § 9 On sotting 

t -1 ^ ) ~ A (if7^)7, r > 

we infer that J^^ lies above but that all <7,^ lie beloiv /. Xho 
limit function f{r, r) is seen to satisfy the integral equation 

/(r,r) = A(/),,,, .(190) 

and to he between the indicated limits t and r+ 1. 


We set J{r,r) = cf{r,r) 

for the solution of J of Problem III and determino the oonstanb 
c as in § 9 The intensity obviously lies between the limits 

c (t + cos 0)<I (t, f) < c (t +1 + cos 0), ... ,(197) 
for all directions r, whence 


hm 


^=0 


uniformly for all directions u This gives, according to tho 
fundamental relation (4*1), 


. (198) 

The proof of the uniqueness of tho solution follows the same 
lines as m § 10 It is sufficient to repeat the mam part of tho proof, 
i e. to show that a solution J of (190) with the properties 


J^Q lini — 
u urn , \ 


Ty~^CX), 


.(109) 


being a given sequence of numbers tending to infinity, vanishes 
identically 

According to the hypothesis made about y, and according to 
(41), we have 


If f 

= 7 - lydoi* ^ const. I oos^ O' dco' 
47rJ J 
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On comparing tins with (44) and (199) we get F~F^ — 0, The 
incident radiation (at t = 0) being zero we thus have, for t = 0, 

F~F^ = ~ f J(O,?)cos0i?a>=:O, 

which IS possible only if the omeigeiit radiation vanishes Thence, 
and from (99), t== 0, wo infer that the Ergiebigkeit must vanish 
too, q e d. 


Remailcs on the uniqueness When the hypothesis made about 
the law of soatteiing is abandoned, Problem III could have 
several solutions This takes place, for instance, in the easily 
integrable limit case, wlioio half the radiation is scattered back¬ 
wards while the other half goes in the oiiginal diiection. Under 
the hypothesis, however, that 


lim 


I{r>r') 
I (t, r) 


holds uniformly for all directions, i.e that the radiation becomes 
isotropic at great depth, the solution is leadily found to bo unique 
Por we find from (41), (44) and (194) that then 


lim 

Tfi=»00 


/(t.O 




holds uniformly for all directions. As in tlio jiioof of Lemma 2, 
§12, we infer, from the positivity ,of the operator, that J s ^Ff» 


§21. DISCUSS CON OP PROBLEM IV 
IN THE CASE (06) 

The optical thickness r of the slab bomg finite, we consider the 
case (66) of constant radiation I* iiioidont at the inner face r=>r*. 

Por the radiation oonnng from above wo use again (99'), while 
the radiation from below is given by (61), or by the more con¬ 
venient formula 

r(T*->T)soofl 

-hi e“®t/(r-hsoosfi,r)ds, 0<7r/2.(200) 
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Theorem XVI Pioblem TV, with the boundary condition 
(66), has a unK^ue solution The Ergiebigkeit lies between the 
hmita _ 1 

.( 201 ) 

the hmits thus being independent of the law of scatteiing 

Pi oof This theorem has, in the ease of uniform scatteiing, 
been proved by Sohwarzschild 

The integial equation (63) of Pioblom IV can be written in 
J (t, ? ) = Ir (J').r.J + ^ )> 

H(t,))=^J . (202) 

= .(203) 


the form 
where 

We set 


In order to compute the quantities 

l)T,r + (t, 7 ), L («^)T,r + K (t, 1 ), 

we first find, fiom (99') and (200), the intensities and I oorie- 
spondmg to the Ergiebigkeiten J i and and then insert them 
into the light-hand side of (41) Aftei some simple computations 
we get 

T -I- cos 0 -H (1 — cos d) e-(T*“T)Boo 0 ^ p <• ,^^2, 
tH-cosO-H I OOS01 e-rl«co0i^ d>7r/2, 


I 

and 


gVJU 


I (t,?) = 


/* 


T*+ 1 


T H-1 -H cos d - COS d e-(T*-T)8oc0^ q 

T-f 1+cosd—(l + oosd)e“’‘i®‘’®^*, d >7r/2. 


These lelations imply, for all directions 7’, the important in¬ 


equalities 


/i(t,7')^Z*--;^^- 


I(t,7')^J* 


T" + 1 

t-H 1-hcosd 
T* -h 1 


On inserting the intensities into the right-hand integral of (41), 
wo obtain, according to 

jeos d'y (t, )•', r) (W = 0, 
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the final inequalities 

L (J,),., + II (r, r) > I* ^ (t, ,).(204) 

L(/),,,+ff(T,0<i*^j = j(T,r) . . ..(206) 

According to the positivity of the operator Ly we see as before 
that tlieio IS a solution between the indicated limits On setting 

J »i+i {ry'i)==^L{J (r, Oj .(206) 

J ^ being defined by (203), we can write (204) in the form 
This imjilies L{J^> L{J-^) and, aocoiding to (206), ?i=l, 2, 
t/g > Jg, and generally ll'ui’thoiinoie, the relation 

^ n) 

and the inequality J—Ji>0 show that generally <7,^< J The 
limit function J of the functions thus exists and has all the 
properties indicated in Theorem XVI. J is, moieover, seen to 
be the V-solution of (03) 

In the present case of a finite t*, the uniqueness is inferred by 
a classical oonoliision, Fust wo note that 


in particular, that 


Q-T\HaoO'\y^l^f „ j 

Itt j ^ 

.(207) 




holds with a suitable constant A. The difforonoe D (t, r) of two 
solutions of (63) satisfies the homogeneous equation D = L{jD), 
Denoting by M the maximum value of [ D |, we find 

|i)|gL(|D|)^ilfL(l)^Ailf; 

in particular M g Ai/, which, according to A < 1, gives M-0, 


Large optical thichneas. In the case of uniform scattering, 
Sohwarzsohild recognized, that for t*->oo, model XV goes over 
into model III. This holds, more generally, for any given law of 
scattering, of the types (11') and (11"). 


IIM r 


5 
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Theorem XVII. For c being fixed, the Ergiebigkeit 

J (T,r) tends towaids cf{r, 1 ) as the optical tluokness t* inoi eases 
indefinitely 

The main jiait of the proof is to show that the Ergiebigkeit 
«7 (r, ?) tends towards a limit function as t*—> co We may content 
ourselves with a brief sketch of the proof Fiom the integral 
eq^uation (63) it can be first recognized that the Ergiebigkoit is 
umformly continuous for all t* Arzela’s selection theoiem then 
guaiantees the existence of one or several limit functions Such 
a hmit function necessarily lies between cr and c (t 1 ). From (11) 
and (202) we have 

I* = ct* e--"* e”" 0, 

as T*-^oo On proceeding, in (63), to the limit t* == co, along any 
special sequence of numbers r*, we find that every limit function 
must satisfy the integral equation of Problem III This equation 
having, however, only one solution cr for r largo, the Ergiebigkoit 
can have only one limit function, as t*“»oo, namely c/ir, r). 


§22 NARROW LIMITS FOR TJIE SOLUTION 
, IN SCHWARZSCHILD’S CASJU, y=l 

In the case of uniform scattering the integral equation becomos 
(64), (64'), the second term on the right being (60), On sotting 


J (t) = I*w(t), 
u becomes the solution of 

u{t)=^L (w)^ + Pjj ('t* - r).(208) 

Let V (t) denote the solution of 

t>(T) = Ir(i;)^+p3(r*-r), (200) 

u IS then the Ergiebigkeit in the case 

•^*=1, .(208') 

while V corresponds to the case 

/(t*0) = cos0; 0< TTj 2, (200') 


On introducing, for any function, the general notation 
*b(T)=cI)(T* —r), 
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we readily find that the operator L, given by (04), satisfies 

L(<h)^ = L(0),*_, = i(cI))^. . . (210) 

ITrom (207), y s 1, we have, aoooiding to (49), = 2, 

1 — L ( 1 )^ + 1^2 + 1-02 ,....( 211 ) 

Furthormoio, we find by means of (79), (81), (82), ?i= 1, 

T = iv(i)^ + —^a + -^(^3 —-Sg) .... (212) 

(208), (210), (211) imply 

whence, aocoiding to the uniqueness property proved in the 
preceding sootion, ,..(213) 

(213) IS physically obvious For U is the Ergiebigkeit when the 
radiation incident is zero at r ~ t* and one at r = 0 + w is there¬ 

fore the Eigiebigkoit when the ladiation incident is one at both 
faces, which yields (213) 

From (209), (211), (212), 

V~V~T~T^‘U, .(214) 

In the sequel wo make use of the 


Lemma (r*,0 )<(t*, 0), 0< 7r/2, implies (r) for 

the 001 responding Eigiebiglcoiten. 

This follows immediately fiom the fact that the Ergiebigkeit 
IS the iY-solution of (64') 

Narrow limits for % (t) can be found by comparison with the 
solution of Problem la. On setting 


'00 

e-«(?(T*-H5)cZs, .(216) 

Jo 


wo infer that and that q^ and q^ increase with incieasing r*, 
having both the limit q^a as r* oo In the model I a, the intensity 
of the radiation coming from below is given by (100) Let us in 
the case F = f denote this intensity by 1, Since q (r) incieases, we 
then find fiom (100) 

T* 4- -I- cos 0<I (t*, 0) < t* + <?2 -i- cos B, .(216) 

for 0 < 7r/2 Taking I as the radiation incident at the inner face 
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r = T* of model la, we fiist see tliat the coiiespondmg Ergiebig- 
keit IS piecisely/(T) = T + g(T) (according to the theoiy of model 
la). On the other hand, the Ergiebigkeit turns out to be 
{r* -f <7i) u + v lesp (t* + q^) u + v, 
when the radiation incident at r = r* is given by the left respective 
light-hand teim in (216) The lemma shows therefoie that, foi 
all T between 0 and t*, 

(t* + <?i) + ^^ </< (t* + (/g) ^ + V, 

which, for the function 

f~f^2r-r*^q-q, 
implies the inequalities 

/-/> {r* + q{}'ii + v- {t* + qz)u-v 
and /"/< (t* + ffg) “ + y <li) ^ ~ 


When, here, ti> and v~v are expressed in terms of u by means of 
(213) and (214), -w is found to he between the limits 


TH-g — g + gi 
T* + ffi + £Z2 
By means of the notation 


< u (t) < 




.(217) 


= .(216') 

Ji 2i 


we can give (217) the foim 

11 \ + + d (r*) 

T* -f 2?iz, (r*) r* + 2m (t*) * 

.(217') 

The use of (217) requires, of course, a numerical knowledge of the 
function q (t) of Problem la. The limits given in (217) are much 
nari'ower than Schwarzsolnld’s limits indicated in Theorem XVI, 
since (/(r) is known to lie between the narrow limits 0-677 and 
0 71 The factor ^(t'*) on the right of (217') is about 0-06 for 
7-* = 0 and decreases very rapidly as r*~^co, 

If, in (217), q is replaced by Milne’s approximate value f, wo 
obtain Milne’s approximate form of J, 




r+f 
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while foi T* large, the tme approximate expression is 

In the case of uniform scatteiing, Theorem XVIT is an evident 
consequonoe of (217) or (217') 

§ 23. THE NET FLUX IN 
SCIIWARZSCIIILD’S PROBLEM 

The flux constant F is not given in Schwarzs child’s pioblem For 
astrophysioal purposes it is, howevei, of importance to know 
Fjl* as a function of t"* We shall now include this function 
between narrow limits 

Fust we must find the analogue of the equation (S4) in 
SohwarzsohikVs model The fiist term on the light must, of course, 
bo omitted, since the radiation incident at the outei face is zeio 
Instead of this term theie is the analogous terra for the radiation 
incident at the inner face, 

Q(T)-=ir oos0e-(T*-^)s‘^o^^J(T* r)c2ai .(218) 

+ 

Furtherinoie, the upper limit of integration in (64) is to be 
replaced by t*. Wo thus obtain Milne’s fiux equation 

Jt Jo 

.( 210 ) 


For T —0, T*, we have 

^F— JF2^dt+ \Q (0), ,,(220) 

Jo 

\F^ - r*JjS?2d/ + ^<3(T*), (220') 

* ® 

whence F= f [J .(220") 

J 0 


Let f resp F' denote the flux oonatant in the case (208') resp. 
(209') Taking account of (49), ?i=3, 4, we note the following 


combinations iii those two cases. 

(3=.2J3g(r*-T), F^f .(221) 

and J^-y, F^F' .(222) 

I 
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In the present case (66) we have 
for the net flux. 

Let us, first, find the relation between F' and P On applying 
(220") to the case (222) we find, accoidmg to (214), 

~ f tF^dt — r* f + (t*). 

Jo Jo 

The flist integral is readily computed fiom (80), = 2, while the 

second one can, by means of (220), taken in tho case of (221), bo 
expressed in terms of P Altogether, the sim^ile relation 

• -(sas) 

IS obtained. 


Tho following remark will be of use The inequality 
i'l(T*,0)<72(T* 0), 0<77/2, 


implies the inequality Fi < F^ for eoiiesponding flux constants 
This follows from tho lemma of tho pieceding section and from 
(218), (220) 

We now lecall the inequality (216), where the incident radia¬ 
tion I gives precisely tho net flux F--^ Taking account of the 
above remark and of (221), (222), we find therefore 

(t*+ ji) J*+J’'< 1 < (t*+(?,)#+ r 


If this is combined with (223), wo obtain simple limits for 

P^Fll\ 4 7 * 

” - * .(224) 


„ 4 7* 


+ "3T*+2g/ 

By means of (216') we can write this 
1 


(r*) 


<ld{T*) .(224') 


This approximate expression of IjP is, even for small t*, very 
aoourato Its value for = 0 turns out to be 0-96, differing only 
6 per cent, fiom the true value IjP ~l. 

In compaiison hereto, we note that the right-hand side of (224') 
gives 8 per cent at t*— 0 For larger r* tho aconiaoy is extra- 
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ordinoiily good According to as wo find, in 

paitioular, from (224'), 

... (226) 

If q is leplaced by its Milne approximation f, wo obtain Milne’s 
appioxiinato formula 

i = 3T*+l. 

The accniaoy is inferred from the fact that 3 w(t*)/ 2 varies 
between the naiiow limits 0 96 and 1 06 



CHAPTER IV 


EXPLICIT SOLUTION OE CERTAIN 
INTEGRAL EQUATIONS 

§24 THE CHARACTERISTIC EQUATION 
We shall now generally tieat the integral equation 

ii{\^~y\)I{y)dy .( 220 ) 

Jo 

by means of Eoimer integrals. The solution will be obtained in 
the form of exphoit integral foimulae. 

Concermng the kernel, we supjiose that, for a certain 5 > 0, 
is quadiatically mtegiablo within 0<a;<oo. Without 
limitation of the generahty we may suppose that 


J [// (a) e®^]2drr < 00; 


s<l 


.(227) 


This hypothesis is certainly fulfilled in Milne’s case. Quadiatio 
mtegiabihty is intioduoed in Older to apiily the l^lanoherel 
Theory of Eourier Integrals 
The solutions of the analogous equation 

r+oo 

/(a;)= H{\x-y\)f{y)dy ....(220') 

J —CO 

are of much simpler form, being aggregates of exponential func¬ 
tions, If w* denote an w-fold root of the ‘ oharaotoristio equation’ 

l = /c(^^)=f'^ Il{\x\)e,^^dt, . (228) 

the funotion (jt) e-~u*v 

IS easily seen to represent a solution of (226'), being an 

arbitraiy polynomial of degree not greater than -1. Wo must, 

of comsG, be sm’e that all integrals involved converge. This is, 

according to (227), obviously true when the real part of u lies 

between the limits i n / , 

— 1 <ic(‘a*) < 1. 

We should expect that the oharaotoristio equation plays also an 
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imporlant role m the theory of (226), and that its solutions show 
the same behaviour for laige x as certain solutions of (220'). 

Wo fiist compute the oharacteristio function. /< {u) in the ini- 
poitant case, wheie II is of the foiin (156), with p (s) increasing 
This includes Milne’s case with p = ^logs k{u) is leadily seen 
to eQ_ual r** r 1 1 ”1 

In Older that (227) be fulfilled, we need only suppose that 
/< (0) = 27/2 (0) be finite. In Milne’s case we have 

. 


the logarithm being zero for -it = 0 
Under the hypothesis (227), the comifiex function t({u) is 
oeitainly holomorpliio in the strip 

|72('m;)1<1. .(231) 

Two obvious pioporties aie 

i< (w) = t< (^), /c (w) = K (— w), .(232) 

showing that k {u) is loal on the imaginary as well as on the real 
axis 

(227) imphes absolute integrability 

J |77(;«)e®*|da3<co; ^< 1, .(233) 

as seen on axiplymg Schwarz’s inequality to (233), the integrand 
being / \ 


(//(a:)e ® '‘)e * 


By simple substitution oj — aj' + TTjt, k {^l) = /c (5 + it) can bo written 
in the form 


whence 


ice 

2|K(i4)|< ■>^(l a; + 2|)_//(|a,|)|e»l»l(2a! 
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If H IS continuous, this inequality shows that 

+ as |i|->c30 .....(234) 

holds unifoimly in every partial strip | s | ^ 6K 1 Theorems of 
Lehesgue, however, show that this still holds under mere measur¬ 
ability oi H As a consequence wo note that the charactoiistic 
equation can have only a fbiito numbei of roots within every 
partial interval | s j g a < 1 

Moie can only be said about the oliaraoteiistio loots when II 
satisfies special hypotheses In the important case, for instance, 
where B is positive, i<{u) is found to increase from /<r(0) to 
><{1) ^ 00 when u goes from 0 to 1 along the real axis (or from 0 to 
— 1) On the imaginary axis k (it) is still real and less than k (0) 
in absolute value The ohaiacteristio equation has therefore the 
double 1 oot rt = 0 when jc (0) = 1 ,. .(235) 

two real roots of opposite sign, however, when 

/<(0)<1</^(1), .(2350 

and at least two imaginary roots when 

/f(0)>l ...(235") 

Besides these it could very well have other complex roots m the 
strip I s I < 1. 

Lotus, now, suppose II to be of the form (156) with increasing p 
In this case, there are, within | Q («) | < 1, no other roots than 
in the case of (236), than two opposite real roots in case of 
(236') and than two opposite imaginary roots in tho case of (236") 
In Milne's case, for instance, all ohaiacteristio roots are thus 
exhausted by the double root w = 0 

/f(w)-/c(w) = 8^s^|^ U=:8 + tl, 

which, according to (232), shows that i< can only bo real on the 
real or imaginary axis Burthormore, tho equation 


I =: 

J 1 




proves that k decreases when u goes from 0 to -l-^co or to — i’co 
along the imagmaiy axis The statement about the roots of /c = 1 
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IS an immediate conseci[uence of these facts It may be added that, 
under the condition (172), i e k (1) = co, the right-hand inequality 
of (135') IS automatically satisfied 


§26. THE INTEGEAL FOEMTJLAE POE 
THE SOLUTION OP (226) 

We shall now construct all solutions of (226) that satisfy the 
condition = ,(236) 

for a) large, a being an arbitral ily fixed number less than one 

T^et "Ug, , ligH 

denote the complete set of chaiacteiistic roots within the strip 
15 I ^ y, If nocessaiy we enlarge y a little such that no loots he on 
the boundaiies | s | = a. We set 

P (lO = (- Uj) (ii - Wg).. {u- Wg,J, ... (237) 

P {u) being oven and real on the real axis, and 


T being holomorphio and free of zoios in the atrip [ s | ^ ot, More¬ 
over, T (it) is even and real for loal w Wo want t (w) in the form 
r (it) “T(. (it)/T„ (it), in such a way that t, (it) is holomoiphic in 
the half-plane — y and that t^{s+U)^ i. as s-^-f-co, while 
T_(5+i0 IS holomorphio for and satisfies t_ ( 5 -i-li)1 as 
s-» — 00 , This obviously umquo repiesontation is, according to 
Cauchy’s formula, explicitly furnished by 


.,(«)=oxp[-^j_-^;;;m*].(239) 


whore ^ > a is chosen such that tho strip [ 51 g /3 contains no new 
roots The logarithm logT(s+iO tho ono that vanishes as 
<->-1-00 Observing that T(it) is never 74ero, and real on tho 
imaginary axis, and that t (« l too) = 1, we soe that t (it) is nooes- 
sanly positive along that axis log t (s -f ii) vanishes therefore loo, 
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when ^->-00 We shall later see tliat the intograls (231)) and 
(239') converge absolutely 
On setting 

o'+('w) = («4+1)”’'^t, (w), o-_(w) = (?t“'1 )”'t_(^0, 

.... ( 2 ‘ 10 ) 

we state the 


Theorem XVIII The integral equation (220) Ims prooisoly n 
linearly independent solutions satisfying (230). 03ioy can bo 
expiessed by means of the forniulao 


^W=<3n-l(«)p'(^|.(2«) 

Qn-i being an arbitrary polynomial of dogroo loss than a, and 

'V27T^Ja-^<a 2 

.(242) 


Corollary The solutions have, for large a;, the form 
fix) ^ SO* (.t) 0 (c'"^”)) 

u* being a oharaoteristio loot within |«|<<y, and IJ* being a 
polynomial of degree less than the mulUpliozty of w*. 

It must be emphasized that, while (220') has piocisoly 2^^ 
solutions with the property (230), (220) has only Jmlf tho mnnbor 
of solutions. 

The abscissa of integration a in (242) can bo moved without 
change of the integral, provided that no polo of r/^ (u) ^ </> [a -1- U), 
1 e no oharaoteristio root is met by it, f{x) vanishes for 0. 


A'pplication to M%lne'‘a caae. Here wo have, according to (30), 

.(2«) 


, . 1 


r^(w) is given hy (239'), while 


^(?i)=:c~r^(^t), .(2.J4) 

c being an arbitrary constant of proportionality. If (242) is 
considered as a Courier integral, wo Mcl by tho Fourier thoorom 

V 27 r Jo 
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This repiesents (j) (it) for all« < 0, since the integral converges and 
IS holomoiphio in the left half-plane We can obtain herefrom an 
exiilicit foimula for the law of darkening in Problem la Prom 
the last foimnla, t=0, and from (188) and (189) we get 


/(O,0) = fJ?’V27rseo0^(~sec5) . . (246) 

The pioportionality constant c m (243) remains to be determined 
such that the coiiespondmg solution/(O/) is for v large Fioin 


(101), I ^0, which gives, together with (246), 


VGtt 


= — lim 3^ (s) 


On the other hand, from (239') and from (244), wo find that this 
limit equals o. We therefore obtain the law of darkening in the 
explicit form 

X (0, 0) = “ jP (1 -1- cos 0) (- sec 6), .(246) 

4 


T_ being determined by moans of (239') and (243). The integral 
ill (239') could be found by numerical integration 

PlancJietd's timrem on XPounei inlegraU. It is convenient to 
state hei e the main theorems of the modern tlieoiy of Fourier 
integrals as far as wo need them m the next sections A complex 
valued function a (<u) of the real voidable a is called quadratioally 
Bummable (q.s.) over —co<a3< +coif 

/» + eo 

|a(a3)pdaj 

j —«> 

is finite. According to Planoheiol, the Fourier transform 

V27rJ-<o 


exists then m the sense of mean convergence, A (f) is again q.s. 
over - c» < ^ < + 00 , and a (a;) is conversely the Fourier transform 

of iT (aj), 1 /»hco 

a (a;)« --p=; A (t) dL 

V2'jtJ -00 
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T^vo q s functions a, b and then Pounoi tiansfornis A, B 
satisfy the Paiseval ielation 

f a{(b)b{x)di~{ A (a) B (.^•) dx 

J —CO J —«> 

In particular, we have 

r+oo r+oo 

\a{x)Y‘dx^\ I -£“1 (x) 12 dx 

J —oo J —00 

The integral (212) is to be understood iii the same souse of moan 
convergence 


§ 26 PUELIWIN'ARIES 

Pi ope} ties of the function r{ic). We know already that r{u) la 
holomoiphic and difPeient from zero within (^xy) 

Moieover, u e know that r (u) convoiges to one as goes to lufnuto 
Within that stiip r (w) being leal on the real and amaginaiy axoa, 
we infer its positivity on these axes. The logarithm of r (-it) con¬ 
sidered in (139) and (139') is therefore real on both axes and toiids 
to zero as u goes to mfimty within the strip logr («) is, furthor- 
moie, an even function of u 


According to (228), k {s-i'it)lV27T is the ITourior transform of 
the qs function a* I), and repieseiits, according to l^lan- 
cherel’s theorem, a q s. function of/, |5 |<1 . (237) and (208) 
show therefore that logr (s + ^/) is a q s function of / too. Wo oaii, 
from this fact, infer the absolute convergence of the Integra la 
(239) and (239') By means of Schwarz's inequality we hud 


47r21 logr^ (w) ja^ | logr (?i) 




I av I 


J/3-10C.’ - > / I t ' J^-ioo I'Ji — 'yp’ 

This shorn moreover, that logT_ (u), being holomorphio m s < S, 
>s bounded in every partial half-plane sS/3'</9 Similar fhota 
are tiue lor log (u) within the half-plaiio s > - 

the“faes™ "" 


I or+(w) I <if 




2 ' 


.(2i7) 
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, , , & 8 
m\u\^<\(T_{u) I < . , (247') 

with suitable xiositive constants m, M 
We mention, finally, that the equation t = t^/t_ can, accoiding 
to (238) and (240), be written 




,(248) 


§27. PROOF OF THEOREM XVIII 
We wiito (240) in the form 

I' +<0 

g{x)^f{x)~ iH\^-y\)f{y)dij, ..,(249) 

J —CO 

whore f{x) = 0, a; < 0, (/ (a) =« 0, x>Q . (260) 

]?or x<0,g (a;) is understood to be defined by the right-hand side 
of (240), Wo now set 

1 r+co 1 /'+« 

.,..(261) 

in other words, wo consider the Fourier tiansforms 
f^(a-fti), y(s-t-^^) 

(legardod as functions of t) of/(a;) e®'® and g (aj) e®'®. 

Since/ (.-r) is supposed to satisfy (236), we may state the obvious 


Lemma 1 , <;^(w) is holomorphio in the half-plane 5< ~a, and 
bounded in every partial half-plane, for instance, iiis ^ — (a -h j8) /2 
From (249), a;< 0, and from (260), 


Il7(a')l^ f 


i « 


^H\!^-'y\)\\f{y)\dy 


< const 

For any A between & and 1, wo find 


r-i« 
0 


/f(|aJ-i/|)|e«V(Zy, 


> l-co ^+00 

\n{\x--y\)\e‘^vdy< \II{\x~-y\)\e^ydy 

0 J -« 

f'+<a 

J —CO 

(/(a;) = 0(e-^i®')) 


whence, for a; < 0, 
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A being an arbitiaiy numbei less than one Consequently, the 
Fourier transfoim has the following property 


Lemma 2 y {u) is liolomorphic in the half-plane a > — 1 and 
boundeclineveiy partial half-plane, for instance, in s ^ — (of -H /2 
It 13 the chief point of the following considerations that the 
regularity half-planes of (j> and u have a common ordinate, 
s = — {y + j3)/2. Along such a common oidmato, wo have, from 
(249), 

1 r+'o r i-w 

y ( 2 t) = 7/(| a:-^ |)/(y )dy 

V 277 J - 00 J ~ CO 

1 /*+«) f I w 

f{y)dy\ II{\x~y\)e»^dx 

W ZtT J ~co J -oo 

hence y («){! —7f (it)}; s=--(<y-l-^)/2 ...(262) 


It is permitted to mterohange the order of integration along the 
common ordinate, since the double intogial converges absolutely 
for a = — (y H- jS)/2 The identity (218) allows us to write (262) m 


the form 


y{u) ^ ((){u) 
o-,,(w) a_{%) 


P{u), 


(263) 


5 = ~ (a-f jS)/2, Here, we know that the left-hand side is holo- 
inorphio for a ^ — (ly -1- j8)/2, while the right-hand side is holo- 
moiphio within a ^ — (a-i-)3)/2. Both sides thus define an entire 
function of u, P (w) being a polynomial of degree 2?i., wo infer 
from (247) and (247') that this entire function is at most of the 
order of w'*' for | u \ large. According to a well-known theorem of 
the theory of complex functions it is therefore a polynomial 
Q (w) of degree not gi eater than n, whence 




Q oan, according to (247'), not bo of degree n since <l>{8 + il)t 
being the Fourier transform of a q,s function, is again a q.s. 
function of t It is thus proved that (j} (u) has the form (241). 

We must not forgot to prove that (241) and (242) actually 
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lepresent solutions of (228) We now start with (241) and study 
the functions/(a;) defined by (242) (241) and (247') show that 
j){%) = 0{\n\~^) for I u I largo, in partioulai, that ^ (s + ti) is a 
q.s function of t Taking account of Lemma 1 and of Cauchy’s 
theorem, we see that the abscissa of integration in (242) may 
be moved to the left without change of the integial Foi any 
5 g (y H-)3)/2, <l>{s + tt) is, Gonveisely, the Fourier transform of 
fix) Hence and from Parseval’s formula, a=6, H = i?, 

«> /* l 

|/('^) l^e^^^dx— [ \^di. 

J —CO J —00 

The right-hand side being bounded for — (a-|-j8)/2, wo infoi, 
on pioceeding to the limit —oo, that/(u) vanishes for, at 
least, almost all negative v (in the sense of Lobesgue) Otherwise 
the left-hand side would increase indefinitely as 5-^—oo 
Lot now, aot = . .(254) 

1 f s y 4- fi 

and (J (.t) = -y=-, y («) dn , s —- 

... (266) 

According to (247), y (u) is regular and 0 {] w |“^) ins ^ — (y -I- ^)/2 
The abscissa of integration can, in (266), bo placed arbitrarily 
far to the right, whence wo infer that g (a) vanishes for almost all 
positive X, 

It remains to prove that f{x) and g{x), found fiom (242) and 
(266), satisfy (249) In other words, wo have to show that, back¬ 
wards, the relation (262), being satisfied by our ^ and y, implies 
(249) We note that the following functions A, B are the Fourier 
transforms of a, 6, 

a iy) == (0 = </> («+ 


My)-e«(^-v)//(|a~y|), B{l)- 


qUx 


K (^t), 


always on the oidmate 5 =~(a+/9 )/2. The Farseval relation 
(a, 5) — (A, gives then, 5 = — (a -1- )3)/2, 

/*00 1 ^*8 I 

/I (I x~yl)/ (y) dy = <!> (it) k (it) du, 

Jo V27ri-J8~'(«) 

.(260) 


11 M P 
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On combining this with (242) and (266), and on taking account 
of (262), we see that (249) and (260), i e that (226) is satiahod 
The right-hand mtegial in (266) converges absolutely becaitso cjt 
and K are both q s functions of t Both sides in (266) loiirosontj 
therefoie, bounded functions of a Tins implies, according to (226), 





The coefficient of the exponent is, hoie, gioater than a Wo 
lemembei, however, that </ was subjected to the only jostriotioii 
that no characteiistic roots should he on ly. Our formnlao 
would, tlieiefoie, not change if a slightly smaller cf be einployocl 
and if ^ also be taken smaller than the original a > There is thus 
no difficulty in establishing (236). 


§28 PROOF OF THE ASYMPTOTIC 
FORM OF THE SOLUTIONS 

According to (241), (li) can have its poles only among the ohar- 
aoteristic roots It has, howevei, at least (?i-i-l) polos within 
151 < a When, in (242), the abscissa of intogiation 5 =: — (a -1- /3)/2 
IS moved to s = (a'f-^)/2, the change of the integral amounts to 
—V2'7r times the sum of the lesiduos of the integrand, oontainod 
in the stiip [ 51 < (/. Near a oharaoteristio root w* of order tho 
integrand is of the form 

e-*" A («) e„-i (u) ; h («*) ^ 0, 

The residue at ' 2 *=has, therefoie, the form 


Q* (rr) 

being a polynomial of degiee less than h. In tho oaso 
<?h_i(w*) 5>^0 the degree equals precisely /c-i. Altogotlior wo 

where/* (a) = 0 (e“i®i) is a solution of (226^), and whero 


V27r^J 



EXPLICIT SOLUTION OP INTEGRAL EQUATIONS 83 

the abscissa of integration now being s = (a + j3)/2 For a>0, 
(226') can bo written 

/*('r)-A(/*)^= f n{'v+\y\)Piy)dy, 

J —00 

A being the operator from 0 to + oo, given by (226). This imiilies, 
according to/* = C) and accoidnig to (233), 

/*(a)-A(/*),= 0(e-% 

foi ^'> 0. The same holds, according to /= A (/), foi the function 
r(r)-A(7%==0(e-“^). ... .(267) 

On the other hand ,) (^’) e®-®, a = (o' + j3)/2, is q s , being the Fourier 
tiansform of along s = (y-h j3)/2. On applying Schwarz’s 
inequality to A (^ )a., the integrand being written 

[//(|r~7/|)c-«i'][r(y)e*«'], 

we infer, theiofoie, that ? (x) = 0 (e~“®), which completes the 
t>roof of the corollaiy. 

§29. NEW PBOOF OF THEOREM X 

The essential hypothesis, for the validity of Theorem X, was the 
Xiositivity of the kernel, 11 > 0, and (119), which can be written 

k(0) = 27/2(0)=1, .(268) 

i e, = 0 IS a characteristic root, being at least double. According 
to T1 > 0, we have, however, 

(0) = 2 J ” OT (0 = 4//j (0) > 0, .(269) 

which shows that w=0 is precisely double. There are no other leal 
roots and, according to | K{it) | </<(0), no imaginary roots. 
Othei complex roots could, however, well exist. The polynomial 
P {%(.) of (237) can thus be wiitton in the form 

P{u)=u^R{u)8{u), 

B (u) = {u~Ui)..,{u~ 

S {u) = {u + W]^) .. (it+ %i~i) • 


where 
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On making, in (241), the particular oboice 
we find 

and, for suitable c, /('») = '« + (?» -1“ o (1) • • .(261) 

(concerning the notation see Theorem VI), since the roots of 
^Sf (^i) “ 0 all he m the right half-plane s > 0, giving rise to lesidues 
that vanish exponentially as a;->0 At the same time it is scon 
that no other choice for leads to solutions of the form (261) 


From (239'), (240) and (260), 

w</»(it)->o, s-^co, .. . .(202) 

and «->0 ..,...(203) 

On tlie othei hand, the first formula (261) is a consequence of 
(242) and holds, here, for all % with 5 < 0, 

rf{x)e0'^dx, .(264) 

V 27 T Jo 

In order to determine c, we observe that (264) implies 

S(^(fi)->>-^; .(206) 


while (261) and (264) imply 


5 V (s) ^; 5 0, 


.( 200 ) 


By combination of (202), (266) and (208), (200), 

SiO) 


A- 




.(207) 


ff_(0) can be found from (239') and (240), 


.(208) 


The abscissa of integration can, here, be moved to the loft and 
deformed into another path, consisting of the parts (- tco, - ip) 
and (A, ^oo) of the imaginary axis, and of tho right half of the 
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ciiole ] -y I = /) T1 i 6 integrals along the first two pieces cancel each 
othei, D“^logT(y) being an odd function The expiession inside 
of the paienthesis, in (268), equals theiefoio 

TT 

~2 


whence, on proceeding to the limit p = 0, the value 

-^logT(O) 


follows. Piom (238), 

r(0) = (-l)«n 

hence 


,c"(0) 21/4(0) 


2/2 ( 0 )^( 0 )“ 8 ^ 0 ) ’ 

a_[0) = (-l)«S{0)lV%aM, 

and according to (267), 

/o=(-i)’-Vm7(o). 

It is to bo noted that ii — 1 is always oven, siiioo the conjugate of 
ohai actons tic lOot is again such a root. 

Ronarh on </«, Computations which may bo left to the reader 
load to the following value of (?«> in (261), 


1 P (2 !<'{% 

J„ 


(U) 


dt. 


In Milne’s case this can be transformed into 


1 p ( 3 tan </i ) 

Jo |sin®<^ f^“tan^) 

The evaluation gives 0*710, 



CHAPTER V 


OTHER PROBLEMS OP 
RADIATIVE EQUILIBRIUM 

§ 30 PURE ABSORPTION NON GRAY MATERIAL 

In this section, a few principal remarks will bo made ooncorning 
the fundamental integial ecj[uation in the case of imiely absorbing 
material with an absoiption coefficient that vanes with the wave¬ 
length. In the case of thermodynamical eq[uihbrium, 

where is a given function of B, 

B,^B,{B)y B^r B,.{B)(ivy .(269) 

JO 

with the following pioperties, 

^^>0, jB„(0) = 0, J3„->oo, B~>co .(269') 

Only these (Qualitative proper ties will bo used in the sequel In 
order to simplify the considerations wo suppose that the abaoip- 
tion coefficient be of the form 

a^(P)=:/5(i/)a(P), .(270) 

The material is again supQioscd to be stratified in paiallol 
planes. The variable t defined by 

T = (^) da; 

J ^00 

is no more the optical depth though being coiivoniont for the 
mathematical discussion. We assume the T-thicknoss of the slab 
to bo finite, t < t*. Tho true optical depth of a Qioint a*, corre¬ 
sponding to v, is 

J —CO 

writing briefly instead of jS (r). 
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The equation of tiansfer is 

^ .... (271) 

while the equation of radiative equilibrium (19) becomes 


COS^ j _Ti 


♦«) ^eo 1 

Jo Jo ‘-iTr 


ludco. 


..(272) 


The net flux ttF of the radiation of all fiequencies is, of couise, 
constant. When (271) is multiplied with cos 0 and integrated with 
respect to ? and v, we find 

K ~ Ft + const., 


K 


JI 


Kydy, ^^V=-0 
0 '^Pv 


IyQo^^0do)i .. . (273) 


which represents a geneialization of Eddington*s relation (32), 
K having, however, no immediate physical moaning in this case. 

Problem V Pind (t, 0) and By (r) from (209), (271) and (272) 
when the radiation incident at the surface t = 0 is zero, while the 
radiation incident at the inner face r — r* has a prescribed value 
ly'^, independentf of the direction (0 < 7 r/ 2 ) 

Wo shall first find the integral equation for B{t), On into- 
giating (271), 

' 1*er^v 800 0{7*~T )+sec q Q-py sooOh-t) ^ | ^ 


ly (r, 0) = 


(3,l600S| e>|. .(274) 


and on integrating, now, through all directions, 


J_ f 

4:7T^ 




I*. 


We introduce the brief notation 


(276) 


fl) (re, P) = 


.(270) 


and 


'¥{B) = jy,B,{B)dv. .(277) 

I This 18 unossontial for llio sixbsoqxiont oonsidoiationa, 
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When the coefficient of ahsoiption is a known funotio 
and frequency, O and 'F are to he consideied definite fi 
their arguments. From (272) and (275), we obtain, i 
mental integral equation 

'j'{jS( t)} = i f"*® {l i--«I.-5(«)}<««+i<?( t*-T 
Jo 

where 0 (a*) [iy 

J 0 

for the determination of JS (r) Once (278) is solved v 
to determine By= (r)} and, from (274), 


§31. DISCXJSSTON OF TTIE NON LXNIOi 
INTEGRAL EQUATION 


The difficulties oonnected with the Pioblom V, 
already in § 4, compel us to content ouiselvos with soi 
tive lomarks on the solution of (278), Again, it la con 
introduce a biief symbol 

J 0 


foi the non-hneax integral oporatoi within (278), (27i 
then 'I'{.B(t)} = 0(.B),+ J(?(t*-t), 

In the special case of gray mateiial, 0 beoomoi 

operator L of Sohwarzschild, while G (a*) simply bocomi 
Tlie operator 0 has a fundamental iiroperty wlii 
called‘monotonity’ and which corresponds to positi 
linear case Since (276) represents, aooordmg to (2( 
creasing function of J?, 


everywhere We also note that the funotion 'F (J3) in 
increasing funotion of B, 


Theorem XIX. Problem V has a solution. 'J’ho int 
tion (278) can be solved by successive apjiroximatioi'i 
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P) 00 f We first study 0{B)^ m the case of a constant B{t) 
TTor any v, B^{r) is then constant too Wo remember that 0{B) 
IS obtained fioin the fiist term on the light of (275) by inultiphoa- 
tion with and by integration through the speotium This teim 


equals, now, r 

wlionoe [“ ^ B^B^ {^,r) dv 

- ijy, BJd^,{T*-T)]dv, ... ( 281 ) 


in the case of a constant B, yVo might, together with (276), 
<I)=<I>j, introduce the functions 


’CO 

<I)„(a.,B)= . (270') 

J 0 


(281) then takes the form 

'F(13) = 0(jS),. + i02(T,i?) + ^(h2(T-^-T,5) .. . .(281') 

ITrom (260') we obtain the following property of the functions 
in particular of Og, 

OgGa, B^co 

This property shows that, for all large enough constants B, 

Oa (t, B)>Q (x), 

0<x< t*, or, according to (281'), 5that the inequality 

^r{B)>0{B)r + lG{r*-T), _.(282) 

0 < T < r*, holds for all sufficiently laige constants B *1 his proves 
the convergenoo of the successive approximations in (278). 
Wo sot xY ^ (^)} = 0 {B,X + (r* - r), .... (283) 

Bq (r) s 0 'F {B) being an inoieasing function of 5, (283) uniquely 
defines a series of functions J5 ,^(t), 7 t== 0 , 1, 2, . Wo obviously 

have A(t)>^oW* 

Supposing the inequality 5,^(r) > JS„_i(t) to be true for one 
particular “W-, we find from (283), according to the monotonity of 
the operator 0, 

'F (t)} > 0 W (^n {'r))> 
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■whence (t) > (r) Induction shows thus the general 

validity of this inequality On the other hand, from (r) < ^, 

Y{5i(T)}<0(5), + iG<(r*-r), 

and accoiding to (282^ -Bi (r) < B. Continuing in this way we 
geneially find (t) < B, The limit function 
JB(T) = limi?,,(T) 

IS as before seen to satisfy the integral equation (278), q gaI. 
The proof of the uniqueness may be omitted, 


§32 REMARKS ON THE CASE OF 
INFINITE OPTICAL DEPTH 


Theoiem I is piobably generally true We shall indicate tho proof 
m the present case 

According to 1^ ^ 0, we find as previously 


,(r,r) = e^‘'se«^i^(r) + jS„sec(?J B^,dty . ^28*1) 


for the radiation from below, 0<7rj2 According to 

0<7r/2, 

and according to (273), 

J v (Ooos^odfco, .(285) 

thus yieldmg 

K{t)> f j{i\T)Qos^Odwy j = ~ 

J + V J 0 p,, 


This IS readily seen to imply i,{i)=.0 for almost all (v, i), booauso, 
otherwise, K (t) would increase exponentially with increasing t, 
in contradiction to (273). * 


Problem VI Pmd B and the inoidont ladiatiou being 
zeio, and the fiux constant J’ ^ 0 being given. 

The integral equation of this problem becomes, according to 
the notation introduced in § 31, 

'F {B (t)}= i f “<I> {| t - « |, jS (t)} dt. 

V 0 


(280) 




OTHER PROBLEMS OE RADIATIVE EQUILIBRIUM 91 

It IS natuial to conjeotiu’o that this equation has a one-parameter 
family of solutions B (t, c) In the simplest case of gray mateiial, 
^^==1, (286) becomes Milne’s linear equation. In the present 
geneial case, the connection between B and F is no longer linear. 
In analogy to the linear case, we could also insert (274), t*=oo, 
directly into the flux integral, whence, according to (276'), 

fOO /'T 

F = 2\ ^l\{t~T,B{t)]dt~2\(l\{r~t,B{l)}di . ...(287) 

When this equation is differentiated and account is taken of the 
relations 

^<V,(a,,J3)= -<!>„(»,£). W{B)=%{0,B), 

(286) is again obtained. On inserting (274), t* = co, into (273), we 
find the relation 

5(i))di=^T + G0nst, 

thus yiotctog 2 

T=00 "I” J 0 

This lelation might be used m order to determine the paiameter 
in B (t, c), when F is given Oiffeientiation of this relation leads 
again to (287). 

§33. ABSORPTION AND SCATTERTNO 
SOIIWARZSOIIILD’S INTEGRAL EQUATION 

When absorption and scattering are simultaneously to be taken 
into account, the determination of the radiation field becomes an 
exceedingly difficult problem It was pointed out in § 4 that tins 
general problem leads to a oomphoated non-linear integral equa¬ 
tion of B, In order to derive this equation we mustfiist, by means 
of (16) and (IC'), express the Brgiebigkeit in terms of the 
emission . This present section is devoted to this partial 

problem. 

Wo suppose scattering to be uniform, ys 1. Since the above 
problem concerns only a definite frequency v, we rejilaoe, for the 
sake of clearness, the affix p by a dash. By,~B\ tf',. . 
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We introduce the total optical depth 

r= I p (y' + a')dx 

J —CO 

and set 


A = ~—,, 

a +a a +ct 


In the outeimost layers of the sun scattering plays the chiof lole, 
while in the deeper layers absorption predominates, 

A~^lj T— ^0, A“^0, T~^oo < <> ^ 288 ) 

The incident radiation being zeio, we obtain, on integrating 
(16) and on inserting the intensity into (10'), 
Schwarzsohild’s integral equation 

J'(T) = A(T)A(J')^ + {1-A(r)}-B'(T) ,.(289) 

for the deteimination of the Brgiebigkeit, A being Milne’s 
oiDeratoi (62) Milne has shown that an equation of the type of 
(289), together with (288), holds under muoli inoie general con¬ 
ditions than under local thermodynamical equilibrium, AA^e con¬ 
tent ourselves with a few qualitative remaiks on the solution. 
In order that (289) have a positive solution, the iV-solution must 
be limto. We piovo that, m general, the Neumann senes cou- 
veiges 


Thbobem XX. Under the hypothesis (288), the W-sories of 
(289) conveiges if ^ q 

holds for some </ < 1. 


P) oof Fust we study the homogeneous equation ooriespouding 
to (289), but with a constant A = A Since (230) lepresents the 
ohaiaoteristio function of Milne’s equation/^ A (/), 



1-m 


will be the oharaoteristio function of 

/(r) = AA(/)^. .(290) 

We know that, tor A = 1, w=0 is the only oharaoteristio root in 
the strip |«| < 1, and that, for A < 1, there are preoisely two roots 
m that strip, being real and of opposite sign. When u increases 
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from 0 to 1, /c (‘It) inoieases from X to + co For an arbitrarily 
given positive j3< 1, we can, therefore, find a 1 such that ^ 
becomes a chaiaoteristio root According to the theory of Chapter 
IV, theie exists a solution of (290), satisfying 

/(T)~e^’-, . .. (291) 

for r large. The following choice of ^ is convenient for our pur- 

(y<j3<l, .. .(291') 

(/ referring to the hypothesis of the iheoiom. 

From (290), 

/+c-AA(/+c)=(l-|)/+c{]-AA(l)}, .. (292) 

c being a constant We observe that A = A(t) is always less than 
one and tends to zero as t->oo Accoiding to the hypothesis 
made about JS' (t), and according to (291) and (291'), the first 
term on the right of (292) will, for all large t, be greater than 
{1 — A (t)} J5' (r) On the other hand we can, according to A (1) < 1, 
choose c so largo that, in the lemammg finite r-raiige, the right- 
hand side of (292) becomes greater than (1 — A) B' This choice of 
c leads thus, for all t, to 

J(T)>A(r)A(/),+{l-A(T)}^'(T), , ...(293) 

J c, J is, of course, positive when c is sufficiently largo, 

This proves, m the usual way, the conveigence of the A'-series. 
Denoting the partial sums of that series by JJ^ (t), e7o = 0, we have 

^ 2i) + (1 ” A) B' 

On subtracting this from (293), 

J-/«,i=AA(J-J„), 

The positivity of A implies thus «/^ < J for all n and for all t. 
From Jq<Ji<Jz<..s, the existence of the limit function J', 
solving (289), IS inferred, q.e d. 

It may be mentioned that the homogeneous equation resulting 
from (289) has no solution (except zeio) piovided that A(t)-» 0 
converges sufficiently rapidly as t-»co. 
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§34 MILiSTE’S MODEL OF A PLANETARY NEBULA 
IN RADIATIVE EQUILIBRIUM 

There is another problem of radiativo oq^uilibinim worth mon- 
tiomng because of diffeient boundaiy conditions A spherical 
gaseous shell is illuminated by a source of light located in the 
oentie of the shell What is the distribution of light in the shell 
when 111 radiative equilibrium? We oonfino our attention to the 
case of puiely absoibing gray material, or to the formally equi¬ 
valent case of monochromatic radiative equilibrium ]jot, for 
the sake of simplicity, the shell be infinitesimally thin, the 
optical thickness t* being kept finite (see tlie beginning of § 6), 

1 e. we neglect ourvatme. 

We introduce the same variables r and 0 as in the previous 
models of a stellar atmospheie The boundary conditions are 
then as follows. 

{a) Theie is no ladiation incident on the outer face, t = 0, 

I(O,0) = O, d>7r/2 

{h) The inner face receives the normally incident parallel 
radiation of the point source 0, the net flux, at t=t*, being ttS, 

(c) This, however, does not ex¬ 
haust the radiation incident at the 
inner face It also receives radia¬ 
tion emergent from other points 
of the inner face, The radiation 
not being weakened while travel¬ 
ling through empty space, the 
figure shows that we have 

1 (t*, tt — 0), d < 7 r }2 

at the inner face. 

This apparently represents a new type of boundary condition. i 
It might be regarded as a special case of 

I (t*, 9) = I (r"^, 7T — 0)'ht{0); 9< .(294) 

% (0), the excess of the incident over the emergent radiation, being 
given. When the finite size of the central star in 0 is taken into 
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account, i (0) lopiesenta a finite function The case of parallel 
radiation noimally incident at t = t* must, as m § 17, be treated 
as a limit case ] 

We now derive Milne’s integial equation of the pioblem, 
under the boundaiy conditions (a) and (204). SohwarzschikVs 
fundamental integral equation of the plane model, with the 
boundary condition (a), and with any i adiation incident at t=t*, 
was givon by (04) and (64'), Accoidiiig to (294), the last term m 
(64') splits into the summands 


... (295) 

The radiation m the hist texm comes through the shell fiom its 
outer parts, 

Z(T*,7r-0)-seo£/ f /(^di, 0<7r/2 

J 0 

Inserting this into (296) and using (49), ?i=l, wo obtain the 
integral equation 

Z'(T) = L*(J).r+:r [ ....(296) 

47rJ + 

whole f'’"* 

L*{J)^mi e/(^)[^(|r-^|) + ^(2T*-T-^)]rf^ 

J 0 

(206') 

On proceeding, as in the beginning of § 17, to the limit case of 
noj mally incident radiation i of net flux ttS, we obtain Milne’s 
equation « 

J(T) = L*(J), + |e--* .(297) 


The net flux ttZ’, being constant in the shell, is found from (294), 


jP = - 

TT 


i (v) cos Odo), 


In the present limit case we find F = S 


I* Foi a finito coniial star, (204) dooa 3iot piociaoly ropi’oaonfc the ooncUtion. o£ 
the pioblom, booauso fclio contial sfcai covers pari of tbo slioll In the limifc case, 
the condition la, of oourao, exact It will, however, bo of advantage m the acquol 
to Qonaidoi tho mathoinatioaUy moio goneial condition (204) 
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It might bo mentioned that the integial equation can easily 
be brought into the simpler foim 

f2r* jS 

J (t) being defined in the larger interval (0, 2t*) by reflection at 
r—r* This shows, as before, that the solution is unique and given 
by the convergent iV-series The same is, of course, true in tho 
more general case (296), since this equation admits of an analogous 
transfoimation 

The radiation emeigent from the outer face is given by tho 
usual formula (47), r = 0 It is of interest to know how tho model, 
S (moie generally ^ {6)) being given once for all, behaves as tiro 
shell becomes moie and more opaque (r^-^-co) Milne’s approxi¬ 
mate formulae suggest that the distribution of light, m particular 
of tho emergent radiation, becomes tho same as in model la, 
F—8 The mere fact that the flux the shell, remains tho 

same for all t*, shows already that the shell cannot become dark 
as T*-^oo The central star is, of course, seen through the shell, 
but its light IS dimmed by absor^ition (or soattoi mg), Tho greater 
part of the light of the shell comes, thoiofoio, from the shell 
itself We now give a rigorous proof of the above-montionod fact 

Theorem XXI When the solution of (207) ap¬ 

proaches the function fjS/(T), 

i e the Ergiebigkeit in the case la, F = 8, tho convergence being 
umform in every finite r-interval. 

Proof, It is convenient to split the equation (297) into two 
simultaneous integral equations, 

= .(208) 

ir being Schwaizsohild’s operator (64), and 

7i(r)=:i)*(/0,+ i rg{t)F{2r*-r^t)cU, .(208') 

Jo 

L* bemg the present operator (296'), Wo obviously have 
' J^g + h; .(299) 
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£7, h depend, of course, upon t* too. The physical moaiung of (298) 
IS that g repiesonts the Ergiebigkeit when the ladiation received 
from the other iiortions of the inner face is neglected. Einally, 
(298') IS the equation of the pioblem, where the tiue boundaiy 
condition (294) is taken account of, i (0) being identified with the 
radiation emergent in the former case 

Accoiding to the syminetiy pioperty ( 210 ) of Schwarzsohild’s 
operator, j/ (t) = ^ (t* — r) becomes the solution of 

!7(T) = ito),+|e-. .. .(300) 

If in the iV-senes leiiresenting g all integrals axe extended ux) to 
infinity instead of to t*, all terms obviously incioase, i.e g (t) 
is smaller than the iV-solution of the same equation, but with 
Milne’s operator A The latter solution has, however, in § 17 been 
shown to be bounded g (r) is thus uniformly bounded for all t*. 
Thcieforo, from (208'), 

h(r)<L*(h)^+jE^(r*-r), ... .(301) 

c >0 being a suitable constant. 

Wo now make a crude oomxiarison with model la. The theory 
of model la shows that there 

1 (t*, 0)~ I (t*, tt — 0) > c' cos 0, 

0<7r/2, c' being a suitable xiositive constant Inseiting tins into 
(200), ,/ being identified with the function /(r) of Problem la, 
wo got 

/(T)>/>(/)r-H-2^3(T*-T) ... .(302) 

(301) iinjilics that h (t) lies below the JV-seiics corresponding to 
the last toim in (301) Similaily, /(t) lies above the JV-sorios 
ooiresponding to the last term in (301) Prom < 2 JS /3 we find, 
thoiofore, „ 

thus yielding, altogether, an inequality 

J(T)<^a(T-|-l), .. ..(303) 

a being independent of t*. 


11 M P 


7 
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(303) Will now bo ubed to show the uniform continuity of J (t) 
foi all Let us flist fix an (othciwiso arbitiary) inteival 
O^T^il We choose t*> h and express, according to (296) and 
(297), the difioience 

J{r")~J{T% OSr'Sr'Sl, 

in terms of integrals Aiianging these intogials in a convenient 
way, we find 


\J{r'')~J{r') 





^’'2 


(H-1) 1 ^ (2t* - t" - 0 - (2t* - r -^ £) I cU 






The second term is loadily seen to bo smalloi than 

gW* (t*-1-1)2 


(t"~/); 




The above inequality signifies, therefore, that in an arbitrarily 
given T-interval (0,/t,), tho fimotions J (r) aro uniformly con¬ 
tinuous, for all T* > /()+1. 

According to Ar 2 ;ola’s selection thoorom, any 8e(|UOiico of 
functions J (t, t*), taken for an arbitrary soq nonce of numbors 
T„*-i-oo, possesses a subsoqucnco that oonvorgos uniformly in 
every finite r-interval Let 3 (t) be tho limit function of siioh a 
Bubseciueiice. It is on account of tho uniform inequality (303), 
that, in (207), the integration and limit jiroooss can bo intor- 
changed, thus yielding 

J(t) = A(J)^, ,7^0, 

A being Milne’s operator Tdieorom V shows, thorofoi’o, that 3 (r) 
is a constant multiple of /(t) Tho oonstaut of proportionality is 
easily determined by means of tho fact that tho fiux constant is 
always iS. On applying the auxiliary thoorom of § 13 to (297) and 
to tho relation 

^ l = i*{l),-|-WW + P.(2T*-T), 
wofind 2 f J(t) [^j(T) + i;j(2T*-T)|dJ = /S(l-e-T*). 
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According to (303), we may again pioceed to the limit t*->oo 
under the integral sign, 

2j“j(T)®a(T) = S, 

thus yielding J {r) — ^Sf{r) Since, thus, all converging sub- 
seq^uences of functions J, r*~>-co, have the same limit function, 
we infer that J (t, t*) converges towaids the mentioned limit 
function as t*->co, the convergence being uniform in every finite 
T-interval, q e d 

When the radiation from other portions of the inner face is 
neglected, i.e when {b) is the only boundary condition at t = t*, 
the integral equation is (298), ^ (r) being the Eigiebigkeit It is, 
in this case, physically obvious that the shell becomes com¬ 
pletely dark as t*~»oo, i e that 


^^0, T^-^cx) ... ,(394) 

This can be rigorously proved as follows, Applying the auxiliary 
thooiem of § 13 to (208) and (300), we obtain 


S f’’* 

— u{t)6~^dl,— g{t)E2{r* — t)dt= 

^ Jo Jo 




Knowing from §22 that u<l and that u-^0 as we 

infer that the last integral m this equation tends to zero as 

In (298), the last term tends to zero as t^-^co, As to the 
integral representing L{g), we spht the integrand into the two 

By application of Schwarz’s inequality, 




Jo Jo 


The first factor has been shown to tend to zero as Ac¬ 

cording to the boundedness of the second factor is loss than a 
constant times f ” J 52 (| ^ |) 

Jo -^ 2(0 

(298) shows thus that (304) is true. 


dt < 00, 
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It should be empliasiaed that this bibliography rofors only to 
mathematical contributions to the problems in question, and to 
the original memoirs concerning tho astronomical and geophysical 
applications. A compiehensivo exposition of tho theory of radia¬ 
tive equilibrium and of the approximate solution of tho problems 
is given in Milne (1), 

To §§ 1-4. The theoretical foundations of tho theory of ladi- 
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ative ecluilibrmm are due to Schwarzsoliild ( 1 ), ( 2 ). He made the 
first fundamental applications to the sun’s atmosphere (law of 
darkening on the sun’s disk> origin of the Fraunhofoi lines in the 
spectrum) To the interioi of a star, the theory was first applied 
by Eddington The mathematical jnoblems arising from the 
steilai intenoi are, however, not taken up m this tract. 

It 18 necessary to say a few woids about the customary intro¬ 
duction of the fundamental quantities of the radiation field The 
use of several difieientials ds, d(Ty dex), dm, d), especially in the 
definition of intensity, is rather troublesome Eioin the point of 
view of the mathematician, it lacks both rigour and beauty The 
entire subject ought to be treated anew, and the fundamental 
equations be derived in a rigorous way, the mam tool being the 
general theory of measure (in the sense of Lebesgue-Hadon). 

To § 6 . The purely absorbing gray model was first introduced 
by Sohwaizsohild ( 1 ), while the model of a purely scattering 
atmosplioio is due to Schuster. The radiation field was originally 
leplaced by two antiparallel stiearns (Schustei-Schwarzsohild- 
approxiniation) The Schuster model was iigorously set up and 
tieated in SchwarzsohikTs fundamental memoir ( 2 ) (case of 
uniform scattering). Einally, the iigoious equations for the ab¬ 
sorbing gray atmosphere were appioximatoly solved by Milne, 
who also studied the spectral distribution of the emergent light 
(cf Milno ( 1 )) 

To § 7 The reduction to integral equations of boundary value 
problems of the elementaly theory of ladiation is duo to Hilbert, 
who used it m Ins attempt to prove Kirohhofi’s laws, The funda¬ 
mental integral equation (28D) for the Eigiebigkeit, obtained from 
(16) and (16') m the piano case, is due to Sohwarzsohild (2). 
The integral equations (63) and (64) of Problems I, II have been 
given by Milne (1), (2) Milne established (63) originally in a 
difieient form, wlnoh is obtained easily by partial integration 
of tlio miegial A (J)^ On writing 

J (0 dl, 

Jo 


C(t) = 


....(306) 
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we find 




= r«(r+0rt^>_ri)e-<*+ 


2t 


j: 


(7 (i 4* t) 

w~~ 


e~^dt 


(306) 


The importance of the positivity of tlio kernels has already 
been realized by Hilbert and Scliwaxzsohild 


To § 9, This method was fiist given in the anthor’s paper (1), 
The existence of a solution could, of couise, also bo obtained from 
SohwarzsohilrVs theoiem about the solution of (64'), the last 
term being (66) Sohwarzschild proved that the solution tends 
towards a limit function as Since (of Theoiem XVt) 

for every t* the solution lies below t +1 and above t, it could 
easily be seen that limit process, t*-^oo, and integration can be 
intei changed in Schwarz schikVs integral equation The limit 
function hes thus between t and t +1 and satisfies Milne's 
equation (cf Kostitzin, Hopf (1)) In view of the central 
significance of Model la, however, a direct treatment seems 
preferable The relation c = was independently proved by 
Bronstem (1) and by the author (6) The simpler proof given here 
is the author’s. 


To § 10. The uniqueness was proved by the author (2). 

To §§11, 12. Concerning the more general equation, treated 
there, see also Hopf (3). 

To §14. Formula (164) for the boundary temperature was 
independently proved by Bronstoin (2) and by the author (6). 
The common root of both proofs is the general formula (161). 
This formula is equivalent with the following faot, proved by 
Amharzumian in the case of Milne’s equation. The resolvent 
kernel K (t, t) of the given kernel // (| r-i |) satisfies the rela¬ 
tion 7^(0, i)=/'(i)//p. Concerning Theorem X see the author’s 
paper (8). 

To §§ 16, 16, The increasing of q (r) was imoved by the author 
(7). The equation (176) is due to Bronstem (3). The numerioal 
lesulte show that, at least for small and largo t, Edchngton’s 
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approximation of q'(t) is the most accurate one (c£ JMilne (1), 
Table I). 

To §17 Model Ih, for parallel incident radiation, was con¬ 
sidered by Milno and Eddington. The case = 0 applies, according 
to Milne (2), to the upper air in the earth’s atmosphere, while the 
case i?* > 0 IS realized in the case of close binary stars (reflection 
efiect) For the approximate solution see Milne (1), 

To §21 Concerning Arzela’s selection theoiem, used also in 
§ (36), see Courant-Hilbert, Methodendei Alathematischen Physih, 
Berlin, 


To §§ 21-23. Theorem XVI is, in the case of uniform scat¬ 
tering, due to Schwaizsohild (2), as well as the conveigonce 
towards a limit model, expressed in Theorem XVII. The 
Sohuster-Sohwarzsohild approximate formula is 




Zii 

t"*-}- 1 ’ 


1 e the moan of Sohwaizsohild’s limits Sohwaizsohild computed 
the coriootioii on replacing the integral equation by a system, of 
linear algebraic equations The approximate expression for J (t), 
given in (217'), can numerically bo coinxmied by using a very 
good appioximation of (t), for instance Eddington’s (Milne (1)), 
which probably differs less than O-S per oont. from the true q (t) 
A oomprehonsivo account of the Sohuster-Sohwarzsohild pro¬ 
blem IS given in Milno (1), § II. 

To §§ 21-30. The solution by moans of Fourier-Laplaco in¬ 
tegrals was given by Norborl Wiener and the author. The 
explicit, though voiy complex, formula for the law of darkening, 
resulting from this theory, is of mam iiitorost in astrophysics. It 
must, liowovor, bo montionod that this method, though solving 
Problem la explicitly, has but a limited apiilicability, as oonoorns 
many other problems of radiative equilibrium. The comparison 
methods, resting on the positivity of the kernel, furnish a simj)ler 
and more general access to a physically satisfactory (thougli not 
explicit) solution. 
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The simplei equation (226'), 

B(t) = J B(t)E(\r-t\)dl, 
J —00 


(307) 


can be interpreted as the approximate foim at great depth of 
Milne’s exact equation. The only characteiistic roots being « = 0 
(double), 5 (t) = ar + & is suggested to be the only solutions of the 
type (236). On using (H05), J = B, the above equation becomes 



C{r + t)~0{T-t) 
2i 


er^dt 


, . (307') 


This equation has been studied by Littlewood and by Haidj'- 
and Titchmarsh, An account of their results is given in Milne 
(1) The latter two authors have proved that the only continu¬ 
ously differentiable solutions of typo (236) are given by 

G (t) = + hr -h c. 


The more general equation, obtained from (307') on reiilacing e- 
by an arbitrary function vanishing like an exponential function 
as <-»oo, has been treated by the author (4), 

The proof of Theorem X, § 29, was given by the author (8) 

To §34. About the relation betwoon radiation and tompoia- 
ture, m absence of local thermodynamical equilibrium, see Milne 
(1), § 19 Oiu X is called there 1/1 

To § 3 6 The ‘ planetar y nebula ’ problem was treated in Miln e 
(3) Reference to the former work by Jeans and Gerasimovich 
IS found at the end of this paper. 



LIST OF SOME FORMULAE 
OF PHYSICAL INTEREST 


Pure absoipUoii. Giay matenal in strict radiative and local 
thermodynamical equilibrium. Model la. 
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ttJP being the net flux, 2’^ being tho temperatnio at the optical 
depth T. y'j 
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Tq being the surface, the effective temperature. 


Tlie function 


TT 

=ip( 


tan^ 




V3<I.(«) = (l + i)oxpj_?J 


sin^ (f) 0 — tan 
21 , log [(1 — GO tan ^)/sin2 


d<l) 


0 u‘^oo8^(}>~\-mn^<j} 

can once for all be tabulated, The law of darkening is (see end of 
§2®) /(O,^)-fF<I)(sec0) 

This holds also in tho absence of local thermodynamical equili¬ 
brium 

Model Ib • I^arallel incident radiation of noimal not flux ttiS, 
tho angle of incidence being O' Tho absorption coefficient for the 
inoidont radiation is supposed to bo a constant fraction n of tho 
general absorption coefficient The following formula© give tho 


upper and inner limiting tomperatiu'o , 2*o 




and 
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Bq = —j- {n SCO 0')\ Ha, = f/S cos Okp pj,sec O' ) 

HJHoo O'fVB. 
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